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Abstract 


-The  objective  of  this  research  was  to  establish  whether  existing 
theoretical  methods  could  be  used  to  predict  the  amount  of  displacement 
that  would  occur  around  8ui  underground  opening  during  or  soon  after  the 
excavation  process.  The  first  step  of  the  study  selected  a  finite  element 
method  of  analysis  based  upon  an  elastic  continuum  as  the  most  promising 
theoretical  method.  The  accuracy  of  the  predictions  was  analyzed  by  com¬ 
parison  with  in-situ  measurements  of  strain  around  as-built  openings. 

For  three  of  the  five  openings  analyzed,  three  appeared  to  present  no  use¬ 
ful  correlation  between  the  magnitudes  of  the  predicted  and  measured  dis¬ 
placements;  for  two  of  the  openings  (Morrow  Point  and  Green  River),  75J» 
of  the  cases  where  the  displacements  were  in  the  same  direction  had  less 
than  an  order  of  magnitude  difference  in  the  amount  of  the  displacement. 

In  general,  a  finite  element  method  based  upon  an  elastic  continuum  does 
not  appear  to  provide  an  acceptable  predictive  method  except  where  the 
rock-system  behavior  approximates  that  of  an  elastic  continuum.  ^ 

/ 
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STRAIN  DISTRIBUTION  AROUND  UNDERGROUND  OPENINGS  - 
COMPARISON  BETWEEN  PREDICTED  AND  MEASURED  DISPLACEMENTS 

1.  INTRODUCTION 

1.1  This  report  summarizes  a  portion  of  the  work  accomplished  on 
Project  No,  DACA-73-68-C-0002.  This  work  was  accomplished  between 

8  September  196?  and  1  January  1971. 

1.2  The  primary  objective  of  this  research  was  to  determine  if  any 
existing  methods  of  analysis  would  predict  the  actual  behavior  of  under¬ 
ground  openings.  A  study  of  the  existing  analytical  methods  Indicated 
that  the  finite  element  method  of  analysis  probably  would  provide  the 
closest  approximation  of  prototyne  behavior.  Therefore,  the  first  step 
was  to  develop  a  computerized  finite  element  method  that  could  use  an 
elastic  continue  as  a  model.  If  the  elastic  continuum  model  did  not 
provide  reasonable  predictions,  then  the  research  would  consider  modi¬ 
fications  that  would  provide  a  closer  simulation  of  the  actuej,  rock 
system. 

1.3  To  evaluate  the  prediction  accuracy  of  the  finite  element 
method,  it  Is  necessary  to  compare  the  predicted  results  with  actual 
measurements  of  displacements  that  occurred  around  actual  underground 
openings.  To  obtain  the  most  valid  such  comparison  and  where  the  vari¬ 
ables  might  be  susceptible  to  analysis,  we  decided  to  study  those  open¬ 
ings  where 

a.  Displacement  measurements  in  the  rock  had  been  obtained  by 
instrumentation  methods  similar  for  all  openings; 

b.  The  more  common  geometries  of  a  circle  and  a  horseshoe  were 
used; 

c.  The  effects  of  different  geological  environments  could  be 
evaluated;  and 

d.  An  appreciable  range  in  the  spans  of  the  openings  was 
represented. 
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1.1*  Using  these  criteria,  the  folloving  openings  vere  selected  for 
study ; 

Straight  Creek  Pilot  Bore 

Horseshoe  (13  ft  span),  Colorado  State  Highway  Dept. 

Pacheco  Water  Tunnel 

Horseshoe  (17  ft  span),  U.S.  Bureau  of  Reclamation 

Navajo  Water  Tunnel 

Circle  (20  ft  10  in  span),  U.S.  Bureau  of  Reclamation 

Green  River  Twin  Highway  Tunnels 

Semi-circular  Horseshoe  (3*»  ft  span),  V^oming  State  High¬ 
way  Department 

Morrow  Point  Underground  Powerplant 

Horseshoe  (57  f t  in  spam),  U.S.  Bureau  of  Reclamation 

1.5  The  initial  stages  of  the  research  also  had  intended  to  study 
one  or  more  of  the  chaonbers  in  the  North  American  Air  Defense  Command 
(norad)  headquarters  in  Colorado  and  the  Veytaux  Underground  Powerplant 
(Switzerland).  However,  examination  of  the  NORAD  data  Indicated  that  any 
compaurative  euialysis  for  this  opening  would  have  to  be  subject  of  later 
research  because  a  three-dimensionad  finite  element  method  of  analysis 
would  have  to  be  formulated;  this  new  method  would  be  required  if  we 
were  to  obtain  a  reasonably  accurate  simulation  of  the  prototype.  As  of 
the  date  of  the  present  report,  such  a  3-D  method  of  analysis  has  not 
been  formulated  (by  us  or  other  researchers)  to  the  degree  where  it  could 
be  used  for  the  intended  purpose.  The  compaurative  analysis  for  the 
Veytaux  openings  had  to  be  abandoned  because  the  designers  amd  construc¬ 
tors  of  this  chamber  were  unable  to  comply  with  their  original  promise 

to  provide  us  with  the  necessary  measurement  data.  The  latter  develop¬ 
ment  was  regrettable  because  the  Veytaux  opening  has  a  spam  of  some  100 
feet,  and  it  is  believed  that  a  partial  failure  occurred  during  construc¬ 
tion.  And,  since  the  opening  was  instrumented  throughout  the  excavation 
period,  the  resulting  data  might  have  been  useful  in  developing  failure 
criteria  for  undergrovmd  openings. 

1.6  In  regard  to  the  latter  situation,  the  prediction  of  safe  spans 
for  underground  openings  still  is  hampered  by  the  lack  of  failure  cri¬ 
teria  for  underground  openings.  To  our  knowledge,  there  has  been  only 
one  tunnel  that  was  properly  Instrumented  at  the  time  of  a  significant 
failure;  unfortunately  these  failure  data  have  been  held  confldentliJ. 
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and  thus  cannot  be  used  for  reporting  purposes.  In  any  case,  such  data 
fr<Mii  only  one  tunnel  would  have  been  insufficient  to  develop  acceptable 
criteria  for  safe  spans.  Because  of  this  lack  of  adequate  field  data, 
oxir  research  has  been  directed  primarily  at  determining  if  the  finite 
element  method  of  analysis  will  give  a  useful  prediction  of  the  displace¬ 
ments  that  will  occur  arotind  underground  openings. 

1.7  After  selection  of  the  tunnel  for  analysis,  the  next  step  was 
to  prepare  all  data  in  a  format  that  would  permit  comparisons  between 
the  analytical  and  the  measurement  plots.  For  this  puirpose  a  computer 
program  was  developed  to  reduce  the  measurement  data  from  the  raw  strain 
measurements  to  a  plot  of  measurement  time  versus  displacements  in  micro- 
inches.  A  major  and  time-consuming  effort  was  required  to  convert  the 
strain  measurements  to  a  format  suitable  for  computer  operation:  although 
all  measxuring  instruments  were  of  the  same  manufacture  and  type,  the 
agencies  doing  the  measuring  used  different  recording  and  collation  meth¬ 
ods.  Thus  many  months  were  required  to  assemble  into  a  common  format,  the 
thousands  of  data  points  that  were  made  available  to  us  on  tape,  punch 
cards,  computer  printout  sheets  and  field  notebooks.  During  this  assem¬ 
bly  of  the  field  data,  finite  element  meshes  were  generated  with  node 
points  located  that  permitted  rapid  extrapolation  of  the  displacements 
expected  at  the  points  where  actual  measurements  were  made  around  the 
underground  openings. 
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2.  THE  INSTRUMENTATION 


2.1  In  the  tunnels  selected  for  study  the  multiple-positicm- 
borehole-extensometer  (MPBX)  vas  used  to  measure  the  displacements  in  the 
rock  around  each  opening  .  It  measures  axial  displacement  components 
within  a  borehole,  and  is  Illustrated  in  Figure  1.  The  instruments  used 
In  the  tunnels  studied  had  a  sensitivity  of  0.000$  in.  and  an  accuracy  of 
O.UOl  in.;  each  instrument  has  3  principal  cosiponents: 

1.  The  Instrument  head; 

2.  The  collar-anchor  assembly; 

3.  The  downhole  anchor  system; 

U.  The  signal  cable; 

3.  The  readout  apparatus. 

2.2  The  MPBX  is  designed  for  installation  in  a  2-l/l(  inch  diameter 
drill  hole  but  can  be  modified  for  use  in  holes  of  larger  diameter.  Gen¬ 
erally,  the  number  of  anchors  used  in  any  one  hole  will  depend  upon  the 
measurement  detail  that  is  desired;  a  frequently  used  system  contains 
eight  downhole  anchors.  Each  anchor  is  spring  loaded  against  the  side  of 
the  hole  and  connected  by  a  0.4l  inch  diameter  wire  to  a  spring-type  can¬ 
tilever  in  the  Instrument  head.  Any  movement  of  the  anchor  relative  to 
the  head  causes  a  comparable  increase  or  decrease  in  the  wire  length  be¬ 
tween  the  euichor  eutd  the  instrument  head;  the  resultant  motion  in  the  csui- 
tilever  is  reflected  by  changes  in  resistsuice  in  the  electrical  strain  gsige 
fixed  to  the  cantilever.  These  changes  are  read  on  a  conventional  readout 
bridge  that  is  attached  to  the  instrument  head  only  during  the  time  that 
the  reading  is  being  made. 

2.3  After  installation,  each  instrument  is  corrected  for  ambient 
temperature  conditions.  The  cantilevers  are  encased  in  a  steel  instrument 
head  to  minimize  possible  damage  from  construction  operations.  Readings  can 
be  taken  periodically,  or  a  continuous  recording  system  can  be  affixed  to 
each  instrument  head. 
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2.U  The  data  we  have  studied  suggest  that  significant  movements  within 
the  rock  system  may  occ\ir  within  a  relatively  short  time  after  installation 
of  the  instrument;  however,  these  movements  may  not  all  be  recorded  when  a 
periodic  method  of  readout  is  used.  This  latter  situation  introduces  the  pos¬ 
sibility  that  a  maximum  strain  may  occur  between  recordings.  Therefore,  a  con¬ 
tinuous  recording  system  is  highly  preferred  during  the  first  week  or  so  after 
installation;  parti ciolarly  during  the  time  period  when  excavation  of  the  work¬ 
ing  face  still  has  a  significant  influence  on  the  strain  at  the  instriunent 
station . 

2.3  Each  MPBX  system  should  be  of  sufficient  length  for  the  last  anchor 
to  be  affixed  in  a  rock  zone  that  has  not,  nor  will  be  influenced  by  the  ex- 

n 

cavation  of  the  opening  .  The  borehole  should  be  so  oriented  that  its  axis 
is  parallel  to  the  anticipated  direction  of  the  movements  most  significant  to 
the  problem  being  studied.  For  this  latter  reason,  when  the  geologic  struc¬ 
ture  of  the  rock  mass  is  such  that  its  behavior  can  be  expected  to  approximate 
a  homogeneous  body,  the  most  desirable  orientation  for  the  instruments  is  to 
provide  one  vertical  MPBX  at  the  crown  of  the  opening  and  a  horizontal  MPBX 
at  each  spring  line  of  the  opening.  In  addition,  for  some  geological  environ¬ 
ments  it  may  be  desirable  to  install  an  MPBX  vertically  downward  to  determine 
the  stredn  adjacent  to  the  invert.  A  major  difficulty  with  this  orientation 
would  be  protecting  the  lnstr\unent  agednst  damage  by  construction  operations. 

If  Information  is  required  on  adjacent  fault  zones,  the  instrument  should  be 
oriented  so  its  longitudinal  ajcis  is  transverse  to  the  principal  planes  of 
movement;  this  permits  measurement  of  separation.  Or  it  should  be  installed 
at  a  small  angle  to  these  planes  of  movement  so  it  will  measure  the  compon¬ 
ents  of  translatory  displacement.  In  a  similar  manner,  the  orientation  of  ma¬ 
jor  joints  should  be  considered  when  locating  the  MPBX's. 

2.6  In  our  project,  52  MPBX's  installed  in  28  tunnel  sections  were 
analyzed.  The  raw  measurement  data  obtained  were  reduced  by  computer  meth¬ 
ods  to  anchor  displacements  relative  to  the  anchor  furthest  from  the  opening 
face;  these  displacements  then  were  depicted  graphically  by  use  of  a  CALCOMP 
plotter.  The  resxiltant  plots  for  the  instruments  studied  are  shown  on  the 
figures  in  Appendix  A.  A  description  of  the  program  used  to  compute  and  pre¬ 
pare  these  plots  is  presented  in  Appendix  B.  Preparing  the  data  in  this 

*In  the  absence  of  specific  data  on  the  influence  of  the  excavation,  present 
practice  is  for  the  final  anchor  to  be  embedded  at  a  distance  equal  to  1.5  to  2 
times  the  opening  diameter.  The  present  study  indicates  the  preferred  dis¬ 
tance  is  2  to  2.5  times  the  opening  diameter. 
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manner  permits  the  quick  recognition  of  what  may  be  anomalous  readings  caused 
by  instrument  faults  rather  than  actual  movement  in  the  rock.  For  example, 
on  Figure  A-19  (2201)  the  large  spike  in  the  center  of  the  plot  for  anchor  5 
probably  indicates  anchor  slippage  or  an  accidental  Jarring  of  the  instru¬ 
ment  and  therefore  shoald  be  disregarded  in  analyses.  The  time- 
displacement  plots  permitted  visual  identification  of: 

a.  the  maximum  amount  of  movement  that  occurred  at  each  instru¬ 
ment  (which  is  designated  by  the  letter  "M"  on  the  compari¬ 
son  plots ) ,  and 

b.  the  magnitude  of  permanent  displacement  after  the  instru¬ 
ment  had  been  positioned  for  a  relatively  long  time  period 
(which  is  designated  by  the  letter  "S"  on  the  comparison 
plots ) . 

2.7  It  is  believed  that  the  "M"  values  will  be  of  use  in  the  design 
of  temporary  supports  for  openings;  the  "S"  values  probably  are  indicative 
of  permanent  splacement  and  thus  would  assist  in  the  design  of  permanent 
supports  and  linings  for  tunnels.  The  "M"  vedues  ailso  could  be  used  to 
estimate  the  amoxmt  of  streiin  (or  displacement)  a  rock  system  csm  accept 
without  immediate  collapse  of  the  opening. 

2.8  As  may  be  apparent  in  the  time-displacement  plots,  it  frequently 
was  difficult  to  establish  the  M  and  S  displacement  values.  However,  if 
continuous  recordings  nad  been  used,  the  peak  smd  the  horizontal  asymptote 
for  each  curve  might  have  been  identified  more  readily  and  accurately.  We 
assumed  that  the  "M"  and  "S"  values  we  selected  would  be  significant  repre¬ 
sentations  of  actual  displacement  that  occurred  in  each  opening.  For  some 
instriunents ,  a  stable  condition  never  appears  to  have  been  reached,  euid 
frequently  it  could  not  be  determined  if  this  was  attributable  to  instru¬ 
ment  faults  or  because  the  opening  actually  dia  not  reach  equilibrium  dur¬ 
ing  the  period  of  measurement.  For  example,  the  plots  shown  on  Figure 

A-T  indicate  that  the  displacements  did  not  stabilize  dviring  the  period 
of  measurement. 

Generally,  we  found  that  the  maximum  recorded  displacement  occvirred 
relatively  soon  after  installation  of  the  instrument  and  that  all  later  dis¬ 
placements  were  a  considerably  lesser  amount.  Here  it  is  important  to 
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recognize  that  because  an  instrument  cannot  be  installed  immediately  af* 
ter  the  excavation  of  the  rock  within  an  opening,  a  significant  amount  of 
displacement  may  occur  that  never  is  recorded  by  this  type  of  instru¬ 
mentation*.  Thus,  where  the  depth  of  cover  permits,  it  would  be  very  de¬ 
sirable  to  install  such  instruments  in  a  hole  bored  from  the  ground  sur¬ 
face  and  have  the  bottom  anchor  in  the  hole  in  close  proximity  to  the 
proposed  excavated  perimeter  of  the  opening;  if  such  instruments  then 
were  attached  to  continuous  recorders  it  would  be  possible  to  measure  the 
amount  of  displacement  that  occurred  during  euid  immediately  after  the  ex¬ 
cavating  process.  Part  of  this  latter  procedure  was  followed  at  the 
Pacheco  Tuiuiel  where  the  MPBX's  were  instedled  from  the  ground  surface 
prior  to  excavation.  Continuous  records  were  not  made,  but  the  records 
clearly  show  that  as  the  excavation  approaches  and  passes  the  instru¬ 
ment,  there  is  a  corresponding  change  in  the  displacement  readings—^. 


^Therefore,  the  values  shown  on  the  plots  in  many  cases  might  be  moder¬ 
ately  to  considerably  incresised  if  the  absolute  displacement  had  been 
recorded. 
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3.  THE  ANALYTICAL  METHOD 


Input  Information 

3.1  Concurrent  with  the  reduction  of  the  measured  data,  the  finite 
element  analyses  were  conducted.  The  results  of  these  studies  were  pre¬ 
dicted  displacements  at  each  of  the  instrument  stations  for  each  geometry 
and  geology  selected.  To  avoid  biasing  the  results,  the  field  informa¬ 
tion  made  available  for  the  finite  element  study  was  limited  to 

1.  Young's  Modulus  and  Poisson's  Ratio  of  the  rock  sis  de¬ 
termined  by  laboratory  tests  on  intact  specimens; 

2.  the  opening  geometry;  and 

3.  the  depth  of  cover  above  the  opening. 

3.2  In  this  stage  of  the  project,  other  geological  factors  such  as 
Joint  systems  smd  faults  were  not  introduced  with  the  exception  of  two  of 
the  instrument  stations  at  the  Morrow  Point  Underground  Powerplsuit.  In 
the  latter  csise,  the  chamber  wsis  intersected  by  two  shear  zones.  The  man¬ 
ner  in  which  these  zones  were  considered  in  the  finite  element  analysis  is 
discussed  below. 


General  Comments 

3.3  Analytical  solutions  to  boundary  value  problems  in  which  gravity 
stresses  are  important,  opening  geometry  is  complex,  and  various  strata 
may  occur  (even  if  the  material  is  assumed  to  be  continuous)  are  difficult 
to  generate.  Such  problems  recommend  the  use  of  numerical  solutions  such 
as  the  finite  element  method-— (  The  finite  element  computer  program  used 
for  the  anetlyses  discussed  herein  is  a  modification  of  that  described  by 
Wilson— {  A  detailed  description  of  the  method  we  used  is  presented  in 
Reference  5* 
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Formulation  of  the  Problem 


3.*t  Each  of  the  underground  openings  and  the  surrounding  material 
are  depicted  as  shown  schematically  in  Figure  2.  The  region  surrounding 
the  opening  is  subdivided  into  a  mesh  of  finite  elements.  Meshes  used 
herein  are  squjure  with  a  height  approximately  12  times  a  characteristic 
dimension  of  the  opening.  Thus  the  top  surface  of  the  mesh  does  not  al¬ 
ways  correspond  to  the  ground  surface.  Boundeory  conditions  applied  to  the 
exterior  of  the  mesh  are: 

a.  The  horizontal  top  surface  is  subjected  to  a  normal 
stress  q  of  magnitude 

q  =  Y  (1) 

in  which  Y  is  the  unit  weight  of  the  rock  above  the 
zone  considered  and  is  the  top  of  the  mesh  below 
the  ground  surface. 

b.  The  vertical  boundaries  are  subjected  to  a  horizontal 
normal  stress.  Increasing  linearly  with  depth  of 
magnitude 

■  K(q  +  yx)  (2) 

in  which  a  is  the  horizontal  normal  stress,  q  is  the 

vertical  normsd  stress  at  the  top  of  the  finite  element 
mesh  as  defined  in  Equation  1,  z  is  the  depth  below  the 
top  of  the  finite  element  mesh,  y  is  the  unit  weight  of 
the  rock  within  the  mesh  and  K  is  a  coefficient  related 
to  the  nature  of  the  prestress  inherent  in  the  rock  mass. 

c.  The  horizontal  bottom  surface  is  restrained  from  vertical 
movement,  but  horizonteLL  movement  of  all  node  points  but 
one,  is  permitted.  This  boundary  condition  is  depicted 
schematically  by  the  roller-supports  shown  in  Figure  2. 

The  one  restrained  node  point,  usually  in  the  middle 

(as  shown  on  Figure  2),  prevents  rigid  body  displace¬ 
ment  of  the  entire  mesh. 
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d.  The  materisd  was  assumed  to  exhibit  a  unit  weight,  y, 
which  was  applied  to  all  of  the  elements  within  the 
continuum. 

3.3  Generally,  a  zero-horizontal  displacement  boundary  condition  is 
used  for  an  outer  vertical  boundary  rather  than  the  stress  boundary  con¬ 
dition  description  under  b.  above.  Such  an  approach  is  useful  if  the  co¬ 
efficient  K  is  equ£d  to  that  for  an  isotropic  linear  elastic  material 
under  plane  strain  conditions,  K^. 


in  which  v  is  Poisson's  ratio.  For  the  more  usual  circumstances  where 
the  rock  has  been  prestressed,  K  is  likely  to  be  larger  than  that  indi¬ 
cated  by  Equation  Results  of  elastic  analyses  for  the  openings  in¬ 

vestigated  were  obtained  for  values  of  K  ranging  from  to  K  =  1. 

3.^  In  calculating  stresses  due  to  the  various  applied  and  gravity 
loadings,  when  an  elastic  anedysis  is  being  performed,  it  was  necessary 
to  apply  all  tht;  loads  simultaneously.  However,  in  order  to  calculate 
the  displacements  within  the  medium  due  only  to  the  construction  of  the 
underground  opening  it  wais  necessary  to  carry  out  two  analyses.  The  dis 
placements  resulting  from  the  excavation  were  considered  to  be  the  differ 
ence  between  these  two  cases (See  paragraph  ^.55). 
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COMPARISON  BETWEEN  ANALYTICAL  AND  MEASURED  DISPLACEMENTS 


4.1  The  points  on  the  comparative  displacement  plots  for  the  five 
tunnels  were  ccnnputed  by  the  methods  described  in  the  preceding  sections 
of  this  report  and  in  Appendix  B.  In  studying  the  plots,  particular  at¬ 
tention  must  be  paid  to  the  displacement  scales:  In  several  cases  the  dif¬ 
ferences  were  so  great  that  it  was  necessary  to  plot  the  measured 
displacements  and  the  ccdculated  displacements  on  different  scales.  In 
these  illustrations,  the  curves  referred  to  as  "C"  or  "K"  are  the  dis¬ 
placement  plots  derived  from  the  finite  element  analyses;  the  "M"  and  the 
"S"  plots  en:e  those  based  on  the  measured  maximum  and  the  measured  "stable" 
or  permanent  displacements,  respectively.  The  latter  points  vere  selec¬ 
ted  as  representative  of  what  appecured  to  be  the  ultimate  and  final  dis¬ 
placement  of  the  instrument;  i.e. ,  there  appeared  to  be  little  or  no 
additional  significant  movement  of  the  anchors. 

4.2  When  the  displacement  is  along  the  axis  of  the  MPBX  toward  the 
tunnel  it  is  plotted  as  "Axis  Extension  (■*'}",  and  when  the  anchor  moves 
away  from  the  tunnel  along  the  MPBX  axis,  the  displacement  is  plotted  as 
"Axis  Contraction  (-)".  On  each  illustration  there  is  a  schematic  pre¬ 
sentation  of  the  orientation  of  the  Joints  around  that  particular  section. 
This  latter  Information  will  assist  in  interpreting  anomalies  that  appear 
in  t.  c  displacement  plots. 

4.2  The  illustrations  are  ordered  according  to  the  span  of  the  tun¬ 
nel  sections  being  discussed.  In  the  initial  discussion  of  the  objectives 
of  this  research,  it  was  indicated  that  we  had  hoped  to  determine  the  in¬ 
fluence  of  the  span  width  and  the  opening  shape  on  the  displacements.  How¬ 
ever,  we  found  there  were  too  many  euiomalies  and  too  few  instrument 
points  to  allow  a  reasonable  analysis  to  be  made  of  the  influence  of  open¬ 
ing  span  and  shape. 


Straight  Creek  Pilot  Bore 

4,4  This  tunnel  ves  constructed  by  the  Colorado  State  Highway  De¬ 
partment;  it  is  about  55  miles  west  of  Denver,  Colorado  on  proposed  High- 

k 

way  1-70.  Its  purpose  was  to  obtain  information  on  the  geologic  and  rock 
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mechajiics  conditions  that  would  assist  in  the  design  and  construction  of 
twin-bore  highway  tunnels-,  the  iatter  are  to  be  located  parallel  to  and 
within  a  few  hundred  feet  of  the  longitudinal  axis  of  the  pilot  bore. 

4.5  The  geometry  of  the  pilot  bore  is  that  of  a  horseshoe  approxi¬ 
mately  13  feet  in  height  and  width.  The  tunnel  is  about  8300  feet  long. 

The  maximum  depth  of  cover  is  about  1500  feet  at  a  point  some  2,000  feet 
from  the  west  portal  of  the  tunnel.  The  actual  amount  of  cover  at  each 
instrument  station  discussed  in  this  report  will  be  noted  in  the  respective 
descriptions  for  these  sections. 

Geology 

4.6  The  bedrock  at  the  surface  and  in  the  pilot  bore  consists  of  the 
Pre-Cambrian  Silver  Plume  granite  with  inclusions  of  older  Pre-Cambrian 
metamorphic  rocks  (Idaho  Springs  formation)  formed  from  sedimentary  rocks 
(Figure  3).  About  75/J  of  the  rock  is  granite,  fine-to-mediiam-grained  and 
consists  chiefly  of  plagioclase  feldspar,  microcline,  quartz  and  biotite 
with  accessory  muscovite.  The  foliation  in  the  rock  is  derived  from  the 
microcline  grains  which  usually  have  a  subparallel  orientation.  The  meta¬ 
sedimentary  rocks  include  a  wide  variety  of  biotite-rich  gneisses,  schists, 
and  locally,  quartzite.  All  rock  is  locally  cut  by  diorite  dikes  that 
probably  are  Tertiary  age.  The  tunnel  is  transected  by  the  Loveland  Pass 
fault  zone.  This  zone  consists  of  numerous  faults  and  shear  zones  of  di¬ 
verse  orientation  and  these  are  separated  by  relatively  competent  ribs  of 
granitic  or  metamorphic  rock.  The  faulted  rock  ranges  frcai  claylike  gouge 
to  highly  fragmented  metamorphic  auid  igneous  rocks.  The  average  attitude 
of  the  faults  and  shear  zones  greater  than  1  foot  wide  is  about  N45°E 
with  dips  of  50°  to  70°SE.  The  average  attitude  of  the  foliation,  particU' 
larly  in  the  meteisedimentary  rocks  is  about  N40°E  with  a  dip  of  30°SE. 

4.7  The  joint  and  fracture  surfaces  commonly  are  coated  and  show 
slickensides.  Typically,  the  Joint  is  coated  with  chlorite  and  the  rock 
for  .01  to  0.1  of  a  foot  to  either  side  of  the  Joint  is  partially  altered. 
The  average  dip  of  the  Joints  is  about  50°  but  they  tend  to  have  a  random 
strike  and  dip. 
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Rock  Properties 


U.8  Results  of  physl co-mechanical  tests  reported  hy  Robinson  & 
probably  represent  maximum  values  measured  on  the  more  competent 
rock  and  are  not  considered  truly  representative  of  the  elastic  properties 
of  the  rock  system.  Table  1  presents  average  values  for  the  granite  and 
metamorphlc  rock  types. 


Instrumentation 

^.9  A  total  of  kh  rock-mechanics  instrument  stations  vere  installed 
in  this  tunnel.  Of  these,  21  representing  9  tunnel  sections  vere  se¬ 
lected  for  study;  the  balance  of  the  MPBX-type  instruments  were  not  con¬ 
sidered  becavise  of  difficulty  in  interpreting  the  measured  displacements. 

In  the  latter  regard,  this  tunnel  was  one  of  the  more  difficult  ones  to 
analyze  because  of  apparent  problems  with  the  instrumentation.  This  was 
one  of  the  first  tunnels  where  extensive  use  was  made  of  the  MPBX  system, 
and  thus  debugging  difficulties  apparently  occurred.  However,  we  believe 
that  sufficient  accurate  displacement  measurements  vere  obtained  to  pro¬ 
vide  a  reasonable  basis  for  comparison  with  the  elastic  analysis.  In  ad¬ 
dition  to  MPBX  and  single-position  extensometers ,  the  loads  on  several  of 
the  steel  supports  were  measured  by  electronic  load  cells  placed  between 
the  legs  of  the  sets  and  the  footblocks,  and  in  a  few  cases  in  horizontal 
positions  in  the  crown  of  the  sets  and  between  leg  and  invert  struts.  The 
data  from  the  load  cells  were  not  considered  in  this  research,  but  they  ar» 
described  in  some  detail  by  Abel'^^ 


Analytical  Model 


U.IO  The  geologic  information  and  rock  property  data  suggested  the 
use  of  a  homogeneous* section  for  the  finite-element  analysis.  The  section 
used  for  this  purpose  is  shown  on  Figure  h.  The  elastic  properties  used  in 
the  analysis  were  an  E  of  8,980,000  psi  and  a  v  of  0.2U3.  This  value  of 
E  compares  favorably  with  the  laboratory  results  and  an  average  E  derived 
from  the  field  sonic  measurements. 


*Homogeneo\is"  in  the  sense  that  the  majority  of  the  rock  was  granite;  in  this 
report,  consideration  was  not  given  to  the  Jointing  and  fa\ilts.  The  influ¬ 
ence  of  such  discontinuities  in  structure  is  the  subject  of  a  later  report. 
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Table  1.-^ 


Rock  Property  (Laboratory) 

Porosity,  percent 
Grain  density,  gm/cc 
Dry  bulk  density,  gm/cc 
Saturated  bulk  density,  gm/cc 
Powder-grain  density,  gm/cc 
Longitudinal  velocity,  fps 
Transverse  velocity,  fps 
Young's  modulus,  dynamic,  psi 
Modulus  of  rigidity,  dynamic,  psi 
Poisson's  ratio,  dynamic 

Unconfined  compression  tests: 

Strength,  psi 

Secant  Young's  modulus,  psi 
Poisson's  ratio 
Tensile  strength,  psi 

Triaxial  compression  tests, 

confining  pressure 
Strength,  psi 

Secant  Young's  modulus,  psi 
Poisson's  ratio 

Triaxied  compression  tests, 

confining  pressure 
Strength,  psi 

Secant  Young's  modxilus,  psi 
Poisson's  ratio 

Final  swell  pressure  of  fault  gouge  and 
bore,  psf 


Granite  Metamorphic 


1.13 

1.62 

2.66 

2.75 

2.63 

2.71 

2.6U 

2.72 

2.6702 

2.7679 

16,976 

15,1*90 

9,699  (■ 

8,267  , 

8.2*5  X  10^ 

7.12  X  10| 

3.38  X  10 

2.76  X  10 

0.22*3 

0.283 

25,600  . 

15,800 

7.57  X  lO*" 

7.33  X  10 

0.22 

0.22 

1,290 

770 

1500  psi 

2*0,600  r 

22*.  900 

9.2*0  X  10 

6.75  X  10 

0.19 

0.17 

2*000  psi 
55,800  (■ 

8.66  X  lo” 

0.21 

altered  rock  samples  from  pilot 

625  -  3,270  230  -  9,000 


Field  measurements  made  in  the  pilot  bore  yielded  a  sonic  velocity 
of  13,750  to  20,150  fps  and  an  electrical  resistivity  of  36  to  2,200 
ohm-m  in  the  undisturbed  rock. 


Comparative  Analysis 

4.11  Station  43  85  (Flgta-es  3A  and  B).  (250'  cover):  The  shear 

zone  appears  to  control  the  displacements,  therefore.  It  was  not  reasonable 
to  correlate  the  results  of  an  analysis  based  upon  a  homogeneous  system 
vlth  the  measured  data.  The  two  figures  show  the  large  discrepancies  that 
occurred. 

4.12  Station  49  27  (Figure  6),  (680'  cover):  With  the  exception  of 

MPBX  2027  the  directions  of  the  predicted  movement  generally  agree  with 
those  actually  measured.  However,  the  analytical  results  suggest  that  the 
sprlngllne  MPBX  orientation  of  45°  Is  an  undesirable  orientation  because 
most  of  the  movement  Is  normal  to  the  Instrument  axis.  Also,  there  would 
be  considerable  Influence  by  the  effective  anisotropy  caused  by  Jointing. 

On  this  basis,  calculations  for  MPBX  2027  show  very  small  movement  but  the 
measured  values  approximate  0.2  in.  However,  it  should  be  noted  that  the 
MPBX  Is  parallel  to  the  Important  Joints.  Also,  it  should  be  noted  that 
the  maximum  displacement  generally  is  away  from  the  tunnel  wall.  This  di¬ 
rection  of  movement  may  have  been  caused  by  blast  damage;  if  so,  it  would 
not  be  predicted  by  the  analytical  methods.  In  general,  to  obtain  agree¬ 
ment  between  the  cedculated  and  measured  values,  ve  would  require  a  K  value 
greater  than  1,  that  is,  approximately  2.  For  the  vertical  MPBX,  the  equiv¬ 
alent  modulus  then  wovild  have  to  be  1/10  of  the  intact  modulus  but  for  the 
nearly  horizontal  MPBX's  the  eqvii valent  modulus  needs  to  be  approximately 
1/100  of  the  intact  modulus. 

4.13  Station  59  +  35  (Flgttre  7),  (1300'  cover):  Note  should  be  made 
of  the  major  perturbation  in  the  time-displacement  curves  that  occurred 
about  one  month  after  the  MPBX  was  installed  (Figure  A-7);  tlds  may  be  the 
result  of  instrument  or  human  error.  The  best  K  values,  that  is  those 
which  have  the  best  agreement  with  the  measured  values,  would  be  1.0  for 
MPBX  2015  and  2021  and  0.321  for  MPBX  2020.  The  very  small  movements  on 
the  measured  plots  suggest  that  the  instrument  was  installed  too  late  to 
record  the  majority  of  the  movement;  the  displacement  shown  for  MPBX  2020 
may  reflect  movement  along  a  Joint  set,  possibly  due  to  blasting. 

4.14  Station  65  94  (Figure  8),  (1300'  cover):  The  extreme  irregu¬ 

larity  in  the  measured  movements  makes  it  difficult  to  obtain  a  fair  com¬ 
parison  with  the  analytical  plots.  However,  it  would  appear  that  the 
closest  approximation  is  for  a  K  value  of  0.321.  One  point  of  interest  on 
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the  measured  plot  of  this  station  Is  the  apparent  occurrence  of  a  so- 
called  "tension"  zone  approximately  1/2  diameter  away  from  the  crown  of 
the  tunnel.  (The  latter  is  discussed  in  paracraph  k.22,) 

U.15  Station  69  +  81  (Figure  9)«  (1300'  cover):  The  most  reasonable 
value  for  K  at  this  section  appears  to  be  0.321  with  an  eq\ii valent  modulus 
of  1/k  of  the  actual  modulus.  Although  the  major  joint  sets  are  approxi¬ 
mately  vertical,  they  probably  have  a  limited  influence  because  of  the 
large  amount  of  cover. 

^.16  Station  73  37  (Figure  10),  (1300'  cover):  A  K  value  of 

0.321  appears  to  give  the  best  approximation  to  the  measured  data.  The 
equivalent  modulus  would  have  to  be  about  1/lOth  of  the  actual  modulus  to 
achieve  a  good  comparison  between  the  measured  and  calculated  data. 

^.17  Station  78  17  (Figiire  11),  (1100'  cover):  For  K  =  1.0  there 

appears  to  be  general  agreement  between  measured  and  predicted  displacement 
insofeir  as  direction  of  movement  is  concerned.  However,  the  prediction 
based  on  a  K  of  1.0  predicts  a  displacement  some  ^0  times  less  than  what 
actually  occurred.  It  is  possible  that  the  large  actual  displacement  may 
have  been  caused  by  slippage  of  joint  blocks  in  this  severely  fractured 
rock  system.  The  equivalent  modulus  would  have  to  be  approximately  one- 
third  of  the  actual  modulus. 

4. 18  Station  8l  +  Ul  (Figure  12),  (1100'  cover):  This  section  was 
in  a  severely  sheared  eu:ea.  The  latter  may  account  for  an  actual  dis¬ 
placement  that  was  some  two  orders  of  magnitude  greater  than  what  was 
predicted  for  a  K  value  of  0.321.  Also,  the  orientation  of  the  gouge 
zone  on  the  right  side  of  the  tunnel  would  cause  it  to  intersect  the 
MPBX;  thus  there  could  have  been  meiss  movement  of  a  large  block  of  rock 
over  seven  of  the  eight  anchors. 

I1.I9  Station  83  93  (Figure  13),  (1100'  cover):  The  large  dis¬ 

crepancy  between  the  measured  and  actual  displacements  probably  ceun  be 
explained  by  the  existence  of  the  clay  gouge  zone  that  intersects  the 
upper  piurt  of  the  tunnel;  this  would  strongly  influence  all  three  of 
the  MPBX's  at  this  station. 

1+ . 20  Station  lOU  +  83  (Figure  iH) ,  ( 68O '  cover ) :  The  K  =  1  value 
appeeurs  to  be  the  most  reasonable  for  MPBX's  200U  and  2005;  and  a  K  of 
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0,321  is  best  for  MPBX  2006.  The  eqiii valent  moduliis  could  vary  as  much 
as  1/lOOth  of  the  measured  modulus. 

U.21  Station  llU  *  $3  (Figure  1$).  (250'  cover):  The  considerable 
discrepancies  between  the  measured  and  actual  displacements  at  this  section 
are  not  readily  explainable  but  may  be  due  to  a  strong  influence  of  joints 
and  their  subsequent  movement. 

The  "Tension  Arch" 

k.22  On  several  of  the  instruments  there  is  a  zone  of  maximum  dis¬ 
placement  that  may  represent  a  "tension  arch"  that  has  been  described  by 
some  authors.  Whether  this  is  a  zone  of  true  tension  ceuinot  be  ascertained 
by  the  type  of  measurements  used  at  Straight  Creek.  However,  it  does  appear 
to  be  a  zone  where  there  is  a  maximum  amoimt  of  displacement  towards  the 
tunnel  and,  therefore,  would  tend  to  be  the  decompressed  area  termed  the 
Trompeter  or  Weber  zone  by  many  European  authors.  The  zone  appears  to  de¬ 
velop  a  maximum  amount  of  "tension"  at  anywhere  from  l/6th  to  approximately 
1/2  the  diameter  away  from  the  opening.  Table  2  shows  the  location  of  this 
zone  for  those  MPBX  stations  at  Straight  Creek  where  it  appears  to  occur. 

The  "tension  arch"  phenomenan  warrants  further  study.  The  difficxilty,  how¬ 
ever,  is  to  design  the  physical  or  mathematical  predictive  model  that  would 
adequately  simulate  this  type  of  stress  distribution.  Such  a  model  appar¬ 
ently  would  have  to  incorporate  the  influence  of  joint  orientation  and  fre¬ 
quency,  ambient  stress  conditions,  and  variations  in  deformation  moduli 
between  different  geologic  strata. 
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Table  2 


Station 

MPBX  No. 

Anchor  No. 

Distance  from  Opening 

1+3  + 

85 

2030 

1 

2 

feet 

1*3  + 

85 

2032 

2  &  3 

3  to  6 

feet 

1+3  + 

85 

2028 

3 

10 

feet 

k9  + 

27 

2026 

2 

4 

feet 

ii9  + 

27 

2025 

2  &  3 

4  to  6 

feet 

65  + 

91+ 

2099 

1 

2 

feet 

69  + 

81 

2018 

1 

7 

feet 

73  + 

57 

2014 

1 

2 

feet 

83  + 

93 

2010 

1 

4 

feet 

83  + 

93 

2012 

1 

6 

feet 

111*  + 

53 

2001 

1 

2 

feet 

111*  + 

53 

2002 

2 

4 

feet 

111*  + 

53 

2003 

1  &  2 

8  to  10 

feet 

^.23.  Attempted  Correlation.  The  magnitudes  of  the  calculated  dis¬ 
placements  appeared  to  have  little  or  no  relationship  to  the  magnitudes 
of  measured  displacements.  However,  it  was  decided  to  determine  if  there 
was  a  statistically  definable  relation  between  these  two  variables.  This 
was  done  by  calculating  the  correlation  coefficient,  r,  between  the  mag¬ 
nitudes  of  the  average  and  the  maximum  measured  displacements  and  the  mag¬ 
nitudes  of  the  average  and  maximum  calculated  displacements,  respectively. 
(In  making  the  analysis  the  one  set  of  apparently  extreme  values  for 
MPBX  2008  were  omitted.)  The  average  values  generally  were  obtained  by 
only  a  visual  approximation  from  the  displacement  plots;  the  maximum 
values  were  scaled  from  the  plots. 

Tne  resulting  values  were  an  r  of  0.202  where  the  variables  were  the 
magnitudes  of  the  average  measured  and  the  calculated  displacements;  and 
an  r  of  0.037  where  the  variables  were  the  magnitudes  of  the  meocimim 
measured  euid  calculated  displacements.  Because  these  coefficients  are 
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Nxmerical  Comparison  Between  Measured  and  Calculated  Displacements 
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Note:  Values  in  (  )  were  omitted  in  determining  correlation  coefficient 

***  The  calculated  values  were  >  than  the  measured  values 

**  Values  not  included  because  measured  value  was  (+)  and  calculated  value  was  (-) 
*  Values  not  included  because  measured  value  was  (-)  and  ceLLculated  value  weis  (  +  ) 


considerably  below  the  5^  level  of  significance  for  the  number  of  data 
sets ,  there  appears  to  be  no  useful  relationship  between  the  aforemen¬ 
tioned  veiriables. 

This  analysis  was  performed  only  for  the  Straight  Creek  Tunnel  data. 
There  were  insufficient  data  points  to  make  a  realistic  anedysis  of  other 
tunnels.  Also,  it  was  not  deemed  advisable  to  incorporate  all  of  data 
from  all  of  the  tiuinels  into  one  correlation  analysis  because  there  was  no 
way  to  accurately  weight  the  influence  of  the  different  tunnel  geometries 
and  geological  environments. 


Pacheco  Tunnel 

Tlus  tunnel  is  being  built  by  the  U.S.  Bureau  of  Reclamation 
for  the  purpose  of  treuisporting  water  on  the  San  Luis  Project  near  Los 
Banos,  California.  The  tunnel  is  a  modified  horseshoe  with  maximum  span 
of  17  ft.  and  a  maximum  height  of  l6  ft.  9  in.  The  tunnel  has  a  length  of 
54,^00  feet. 

Geology  (Figure  l6) 

U.25  The  two  major  rock  types  at  the  site  are  Cretaceous  metasand¬ 
stone  and  metashale.  The  major  faults  are  steeply  dipping  and  intersect 
the  tunnel  axis  at  nearly  right  angles.  Shear  zones  associated  with  the 
faults  range  up  to  ?.  ft.  in  thickness.  Groundwater  inflow  during  the  ex¬ 
cavation  ranged  between  2  and  Uo  gpm.  Much  of  the  tunnel  rock  is  severely 
fractured;  some  of  the  rock  contains  quartz  veins  with  a  porous  or  vuggy 
structure.  The  foliation  is  contorted  and  there  is  considerable  variation 
in  the  dip  of  the  Joints.  At  Station  132  +  50  the  metasandstone  strikes 
I'i50°W  and  dips  30°  to  55°NE.  The  bedding  is  thin  to  thick  and  dip  of 
joints  usueLIly  is  greater  than  75°.  The  Joints  are  coated  with  yellow- 
brown  oxides.  At  Station  151  +  68.1  the  metasandstone  strikes  from  N20° 
to  60°W  and  dips  5°  to  20°NP;.  The  bedding  is  thick  to  massive;  the  Joints 
are  spaced  1  to  3  feet  and  have  some  surfaces  that  are  weatliered  grey-brown. 
The  Joints  adjacent  to  fa\ilts  have  a  yellowish-brown  clay  film.  There 
is  10  to  20%  metashale  interteuued  with  tne  sandstone  and  this  contains 
both  calcite  and  quartz  in  veinlets  up  to  a  1/2  inch  in  thickness. 
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Rock  Properties 

k.26  It  vas  difficult  to  obtain  accurate  tests  on  the  rock  because 
of  the  severe  foliation,  the  quartz  stringers,  emd  frequent  fracturing. 
However,  the  results  of  static  and  dynamic  tests  performed  at  the  Purdue 
University  rock  mechanics  laboratory  gave  the  results  shown  in  Table  U. 

Table 

ROCK  PROPERTIES  (Physico-Mechanical)  and  Types  of  Tests  Used  (laboratory); 
Static  and  dynamic  tests.  (Tests  by  rock  mechanics  laboratory,  Purdue 
Ifaiversity. ) 


Specimen 

No. 

Rock  Type  Young's  Modulus, 

E,  X  10^  psi 

Poisson's 

Ratio 

Density, 

pcf 

1 

Metasandstone 

6.2 

0.318 

173 

2 

Metasandstone 

5 

0.21 

3 

Metasandstone 

h.9 

0.39l» 

Instrumentation 

k.27  Two  MPBX's  were  installed  in  boreholes  drilled  vertically  from 
the  surface  prior  to  excavation  of  the  tunnel.  Because  of  severe  fractur¬ 
ing  in  the  rock,  the  drill  hole  containing  MPBX  2208  was  cemented  and  re¬ 
drilled  to  Insure  positive  anchor  setting. 

Analytical  Model 

k.2Q  The  elastic  model  used  in  the  finite  element  analysis  is  shown 
in  Figure  17.  The  rock  properties  used  as  input  were  an  E  of  5*000,000 
psi  and  a  v  of  0.20.  The  depth  of  cover  at  both  sections  was  90  ft. 

Comparative  Aneilysis 

4.29  The  measurements  plotted  on  these  two  sections  (Figures  I8  and 
19)  are  of  considerable  interest  because  they  record  the  displacements  as 
the  excavation  approached  the  instrument.  These  were  the  only  MPBX's  in 
the  entire  study  that  had  been  installed  prior  to  actual  excavation. 
Therefore,  it  is  reasonable  to  assume  that  the  measured  displacements 
represent  absolute  values  of  displacement 
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U.30  The  large  difference  between  the  calculated  and  the  measured  dis¬ 
placements  is  not  unexpected  because  a  tvo-dimensloned  static  analysis  cannot 
reflect  the  behavior  of  this  rock  system  as  the  tunnel  approaches  the  instru¬ 
ment.  Also,  it  is  believed  that  for  Pacheco  rock  the  modulus  of  elasticity 
of  the  intact  rock  tested  in  the  laboratory  is  considerably  higher  than  the 
modulus  of  the  rock  system  in  situ. 

Mavajo  Tunnel  No.  1 

U.3I  This  tunnel  was  constinicted  by  the  U.S.  Bureau  of  Reclamation  on 
the  Navajo  Indian  Irrigation  Project,  San  Jueui  County,  New  Mexico.  The  tim- 
nel  is  10,079  feet  long  and  has  a  circular  cross-section  with  a  diameter  of 
approximately  21  feet.  This  tunnel  was  excavated  full  face  by  a  boring 
machine. 

Geology  (Figtire  20) 

U.32  The  formations  at  the  site  are  alternating  beds  of  Cretaceous 
sandstones  and  shales.  The  structure  generally  is  massive  and  there  is  no 
evident  Joint  system  although  at  Station  I6  +  there  were  randomly  ori¬ 
ented  sliekenslded  Joints  in  the  shale.  The  sandstone  is  fine  to  medium 
grained  with  occasional  pockets  of  medi\un  to  coarse  sandstone.  There  sure  oc¬ 
casional  lenses  of  silty,  slightly  limey,  carbonaceous  or  shaley  sandstones. 
The  rock  generally  is  moderately  hard  to  moderately  friable  and  poorly  lam¬ 
inated  to  cross-bedded.  The  shale  is  silty  to  clayey,  slightly  carbonaceous, 
partly  sandy,  moderately  hard  and  laminated  and  fissile  to  blocky.  There  is 
a  transition  material  of  shaley  sandstone  that  is  fine  to  medium  gredned, 
silty,  carbonaceous  and  moderately  heurd  to  moderately  friable  and  is  lamin¬ 
ated  to  poorly  laminated.  The  beds  generally  have  a  very  lov  dip  or  are 
horizontal.  A  continuum  model  would  appear  to  closely  approximate  this  geo¬ 
logical  environment. 

Rock  Properties 

h.33  Laboratory  tests  gave  secant  E  values  at  a  first  loading  of  1000 
psi  of  507,000;  567,000;  and  572,000  psi  for  the  sandstone.  Its  compressive 

0/ 

strength  averatged  3,010,  5,730  and  1,820  psi,  respect! vely^( 
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Instrumentation 

The  Instruments  consist  of  one  MPBX  at  eacV  of  three  stations  and 
convergence-measxiring  plugs  at  two  stations,  (information  on  the  latter  may 
be  found  in  Reference  2.) 

Analytical  ^fodel 

k,35  The  elastic  finite  element  model  is  shown  in  Figure  21.  Rock  prop- 

erty  input  values  were  an  E  of  500,000  psi,  v  of  0.2  ,  and  a  density  of  l60 

pcf.  The  depth  of  cover  at  all  three  stations  was  350  ft. 

Compeirative  Analysis 

U.36  Station  7  (Figure  22).  The  MPBX  was  installed  in  an  old  head¬ 
ing  that  had  been  excavated  1-1/2  years  prior  to  the  installation  of  the  instru¬ 

ment.  This  long  delay  probably  resulted  in  most  of  the  stredn  being  relieved 
prior  to  installation  of  the  instruments. 

U.37  Station  lU  +  76  (Figure  23).  This  MPBX  was  installed  l6.5  ft  be¬ 
hind  the  advancing  face.  Since  the  tunnel  was  being  driven  full  face  by  a 
boring  machine,  the  instrument  probably  was  installed  within  2U  hours  after 
the  machine  had  passed  the  Instrument  point.  Despite  this  rapid  installation 
of  the  instrument,  the  amount  of  measured  displacement  is  so  small  as  to  ap¬ 
proach  the  accuracy  of  the  instrument,  and  for  that  reason  may  not  represent 
a  true  measure  of  the  absolute  displacement. 

4.38  Station  l8  85  (Figure  24).  This  instrument  was  installed  19.3 
feet  behind  the  advancing  face;  the  latter  means  that  it  probably  was  installed 
within  24  hovirs  of  the  time  the  tunneling  machine  had  passed  the  instrument 
station.  Although  there  is  an  order  of  magnitude  difference  in  magnitude,  the 
directions  of  the  measured  and  the  calculated  displacements  are  the  same.  The 
lack  of  comparison  with  Station  l4  *  76  is  not  readily  explainable. 

4.39  An  overall  interpretation  of  the  instrument  behavior  at  Navajo  is 
that  the  rock  may  have  behaved  as  an  elastic  continuum;  if  so,  there  wovild  be 
an  almost  complete  release  of  strain  prior  to  the  time  the  instruments  were 
installed.  If  this  assumption  is  true,  it  then  supports  our  earlier  state¬ 
ments  that  for  accurate  displacement  measurements,  it  is  essential  for  instru¬ 
ments  to  be  installed  prior  to  the  actual  excavation  of  an  opening;  the  latter 
particularly  would  be  true  when  the  behavior  and  composition  of  a  rock  system 
approaches  or  equals  that  of  an  elastic  continuum. 

*In  the  absence  of  a  laboratory  value  for  v,  this  value  was  asstuned  as  repre¬ 
sentative  of  the  sandstone. 
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Green  River  Tunnels 


4.1*0  The  Green  River  Tunnels  were  excavated  by  the  Wyoming  Highway 
Department;  they  are  twin  tunnels  for  U.S.  Highway  I-80  through  Castle 
Rock  Ridge  near  Green  River,  Wyoming.  Each  tunnel  is  a  horseshoe  shape 
with  an  excavated  diameter  of  approximately  34  feet,  a  height  from  invert 
to  springline  of  9*66  feet,  and  a  length  of  1100  feet. 

Geology  (Figure  23) 

1*.1*1  The  rock  in  the  tunnels  consists  of  marlstone,  silt  stone, 
silty  sandstone  and  shale  of  the  Tertiary  Green  River  formation.  The 
beds  are  predominantly  flat  with  local  1°  to  4°  dips  resulting  from  gen¬ 
eral  folding  on  northerly  trending  axes.  The  maximum  overbmrden  is  ap¬ 
proximately  250  feet.  The  rock  varies  in  hardness  with  the  shales  being 
the  weakest  and  bedding  planes  tend  to  be  weak.  Prolonged  exposure 
causes  some  local  slacking  but  does  not  produce  excessive  disintegration 
or  volume  change.  Generally  the  tunnels  were  dry  althouj^i  there  was  a 
small  amount  of  seepage  near  the  west  end  of  the  east  bound  tunnel.  The 
Joint  sets  have  vertical  dips  and  are  oriented  N22°W  and  I/80°E. 


Instrumentation 

4.43  At  each  of  the  four  instrument  stations  there  were  prop  load 
cells  and  three  MPBX's.  The  MPBX's  had  a  reunge  of  0  to  1.0  in  and  a 
sensitivity  of  .0002  inches. 
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Anadytleal  Model 

U.UU  The  elastic  finite  element  model  used  for  our  analysis  is 
shown  in  Figure  26.  The  rock  properties  used  as  inputs  were  an  E  of 
150,000  psi,  a  V  of  0.1,  and  a  density  of  130  pcf.  For  this  analysis 
the  depths  of  cover  for  the  various  sections  were  rounded  off  to  90  ft 
for  98  +  35,  2U0  ft  for  103  +  51  and  105  +  05  and  208  ft  for  105  +  5**. 

Comparative  Displacements 

U.U5  Station  98  +  35  (Figure  27),  (80'  cover) (South  Tunnel): 

Ihere  is  relatively  good  agreement  in  both  magnitude  emd  direction  of 
displacement  between  the  measured  and  the  calculated  values;  the  excep¬ 
tion  to  this  is  that  the  movements  of  MPBX  2108  anchors  6  to  8  are  in 
the  opposite  direction  from  that  calculated  for  these  smchors.  There 
does  not  appear  to  be  a  reasonable  explanation  for  this  contradiction. 

It  is  of  interest  that  a  "tension  arch"  appears  to  have  developed  be¬ 
tween  anchors  2  and  U  and  at  anchor  6  in  MPBX  2107  and  between  anchors 
3  and  6  and  anchor  7  in  MPBX  2106. 

U.l»6  Station  103  ♦  51  (Figure  28),  {2k'J'  cover) (South  Tunnel): 
There  generally  is  good  agreement  in  the  magnitude  of  the  calculated  and 
the  measured  displacements.  However,  there  is  a  reversal  in  direction 
in  the  measured  data  on  MPBX  2101.  The  calculated  displacement  using  a 
K  of  .11  appears  to  provide  the  closest  simulation  to  the  measured 
points.  MPBX *8  210i»  and  2105  both  indicate  the  possible  development  of 
a  "tension  arch"  at  anchors  3  through  5  in  210U  and  anchors  5  and  6  in 
2105.  The  movement  of  the  tunnel  crown  away  from  the  opening,  as  demon¬ 
strated  by  the  measurements  in  MPBX  2101,  is  not  readily  explainable. 

U.U7  Station  105  05  (Figure  29),  (233'  cover) (North  Tunnel): 

There  is  relatively  good  agreement  in  the  direction  of  displacements  be¬ 
tween  the  meMured  and  calculated  values.  However,  the  magnitudes  vary 
considerably.  It  is  of  particular  interest  to  note  that  for  MPBX  2100 
the  reversed,  in  direction  of  displacement  eis  shown  by  the  calculated 
plot  has,  except  for  anchors  1  and  2,  a  comparable  directional  change  in 
the  measured  displacements. 


U.li8  Station  105  (Figure  30),  (208'  cover )( South  Tunnel) : 

There  is  generally  good  agreement  between  the  orientation  and  the  magni¬ 
tudes  of  the  measured  and  the  calculated  displacements  except  for  MPBX 
202^  where  the  measured  displacement  is  in  an  opposite  direction  from 
that  csdculated. 

^.^9  In  general  there  is  good  agreement  between  the  measured  and 
calculated  displacements  for  the  measurements  at  Green  River  Tunnel. 

The  anomalies  probably  can  be  explained  by  observation  of  the  time- 
displacement  plots  which  appear  to  indicate  large  movements  occurring 
within  a  relatively  few  hours;  therefore,  a  significant  amount  of  move¬ 
ment  probably  occurred  before  the  instruments  had  been  installed.  A 
major  reason  for  the  relatively  good  agreement  between  the  measured  and 
the  calculated  data  is  that  this  tuimel  is  located  in  a  geologlccd  environ¬ 
ment  that  has  had  only  relatively  modest  amounts  of  orogenlc  disturbance 
and  there  are  no  major  faults  or  shear  zones.  Thus,  it  is  reasonable  to 
assume  that  the  in  situ  behavior  would  approximate  that  of  an  elastic 
continuum. 


Morrow  Point  Underground  Powerplant 

^.^0  The  chamber  for  this  plant  was  constructed  by  the  U.S.  Bureau 
of  Reclamation  on  the  Gunnison  River  some  20  miles  east  of  Montrose,  Colo¬ 
rado  and  adjacent  to  the  Morrow  Point  Dam.  The  chamber  is  206  ft  long, 

57.^  ft  wide,  65  to  13^  ft  high  and  is  about  UOO  feet  beneath  the  ground 
surface. 

Geology  (Figure  31) 

U.5  The  chamber  is  in  a  region  of  anticlinal  folding  and  is  intersec¬ 
ted  by  several  shear  zones.  The  chamber  is  Intersected  by  a  major  shear 
zone  that  has  a  strike  of  N28^W,  a  dip  of  3^^NB,  and  an  average  thickness 
of  2.5  feet.  A  second  major  shear  zone  Intersects  the  upper  esust  end  of 
the  chamber  and  lies  above  the  remainder  of  the  chamber.  The  orientation 
of  the  three  major  Joint  sets  that  intersect  the  openings  sure:  N37°E, 

80°NW;  n63°W;  80°SW;  and  Nl*°E,  1»2°SE.  The  rock  types  primarily  are  bio- 
tite  and  mica  schist,  quartzite  Euid  pegmatite. 
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Rock  Properties 

U.32  Laboratory  emd  field  tests  were  performed  to  determine  physico- 
mechanical  properties  of  the  rock.  A  summary  of  the  results  of  these  tests 
are  shown  in  Table  6. 


Table  6 


8.10/ 


ROCK  PROPERTIES  (Physico-Mechanical)  and  Types  of  Tests  Used  (in-situ  and 


laboratory) : 

Rock  Rock  lype 

Unit 

Young's 
Modulus,, 
E,  X  10^ 
psl 

Poisson's 

Ratio 

Ultimate 

Strength, 

psi 

Number 

of 

Tests 

A 

(Biotite  Schist 

0.86 

0.02 

3,615 

2 

A 

(Mica  Schist 

1.20 

O.Oh 

6,385 

2h 

■D 

(Micaceous  Quartzite  U.lU 

0.06 

15.590 

20 

a 

(Quartzite 

8.81t 

O.lU 

28,820 

k 

c 

Pegmatite 

2.68 

0.05 

17,700 

h 

Instnunentation 

U.53  Precise  leveling  techniques  were  used  to  determine  the  absolute 
displacement  of  25  points  on  the  rock  arch.  Movements  of  the  walls  of  the 
excavation  were  measured  for  nine  pairs  of  points.  Permanent  reference 
points  were  established  above  and  below  shear  zone  A  at  two  stations  to 
determine  if  excessive  movement  of  the  A-line  rock  wall  was  confined  to 
the  shear  zones.  Twenty-five  instrumented  rock  bolts  were  installed  (20 
feet  long)  in  the  walls.  7  MPBX's  were  installed;  3  in  the  arch,  2  at 
the  springline,  and  2  in  the  walls.  Two  single-position  borehole  exten- 
someters  were  Installed  to  a  maximum  depth  of  115  feet.  The  displacements 
shown  for  the  sensor  head  ("anchor"  O)  were  obtained  from  the  plots  in 
Reference  11;  all  other  measurement  data  were  supplied  directly  by  the 
U.S.  Bureau  of  Reclamation. 
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Analytical  Model 

I1.55  The  elastic  finite  element  model  is  shown  in  Figure  32.  Three 
types  of  materials  appeared  to  have  significamtly  different  elastic  pa¬ 
rameters;  the  properties  eulopted  for  these  are  as  follows: 

Rock  Type  A:  £  =  2.5  x  10^  psi 

B:  E  =  3.5  X  10^  psi  V  =  0.03 

C:  E  =  0.1  X  10^  psi* 

The  mesh  in  Figure  32  was  used  for  the  initial  calculation  and  displace¬ 
ments  prior  to  the  removal  of  material.  The  second  mesh  in  which  the 
elements  in  the  excavated  area  were  removed  was  used  for  the  second  com¬ 
putation.  The  area  of  particular  interest  in  Figure  32  is  that  zone  in 
which  the  elements  are  smallest. 

Comparative  Analysis 

1*.55  MPBX's  2201  and  2202:  (Figures  33  and  34).  The  considerable 
axial  extension  between  anchors  1  and  6  on  these  instruments  may  be  in¬ 
dicative  of  the  so-called  "tension  arch".  For  MPBX  2201  the  measured 
displacement  is  less  than  the  calculated  displacement;  for  MPBX  2202  the 
calculated  average  displacement  is  less  than  the  measured  average  dis¬ 
placement.  These  contradictions  are  not  readily  explainable.  It  should 
be  recognized,  however,  that  these  instruments  are  in  the  hanging  wall  of 
the  shear  zone  and  do  not  appeeur  to  be  intersected  by  any  fault.  Any 
movement,  therefore,  should  be  due  to  loosening  of  Joints  and  fractures 
in  the  rock.  In  the  one  case  (MPBX  2202),  it  is  suggested  that  a  con¬ 
siderable  part  of  the  elastic  strain  was  not  released  during  the  period 
of  meets urements;  whereas  for  MPBX  2201,  a  considerable  portion  of  elastic 
strain  may  have  been  released  prior  to  installation  of  the  instrument. 
However,  of  most  vedue  is  that,  with  the  exception  of  the  sensor  head, 
the  direction  of  the  measured  displacement  is  in  agreement  with  that  pre¬ 
dicted  by  the  calculated  displacements.  No  reMonable  explanation  could 
be  found  for  the  apparent  inconsistency  in  the  direction  of  movement  of 
the  sensor  heads  for  these  MPBX's. 

•This  "C"  material  represents  the  shear  zone  material  Euid  is  not  the  same 
rock  xmit  as  "C"  in  Table  6. 
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1».56  MPBX  2215:  (Flg\ire  35).  In  this  section  the  measured  movement 
vas  in  the  same  direction  but  several  times  greater  than  that  predicted. 
This  sharp  increase  in  mag;nitude  may  have  been  caused  by  the  shear  zones 
intersecting  the  instrument.  This  latter  conclusion  is  supported  by 
other  field  observations  that  indicated  a  large  wedge  of  rock  in  this 
peurt  of  the  wall  vas  moving  toward  the  chamber,  and  its  boundary  vas  be¬ 
tween  anchors  6  and  If  the  latter  assumption  is  correct,  then  this 

type  of  inelastic  displacement  would  not  have  been  predictable  by  our  use 
of  an  elastic  model  in  the  finite  element  analysis. 

4.57  MPBX's  2203,  2206,  2207  and  2217:  (Figure  36).  The  wide  vari¬ 
ance  between  the  measured  and  the  calculated  displacements  for  MPBX  2203 
may  be  attributed  to  the  same  situation  described  for  Figure  33  and  34; 
i.e.,  loosening  of  Joints  in  the  rock  would  have  caused  considerable  dis¬ 
placement  and  the  influence  of  Joints  vas  not  introduced  in  the  finite 
element  analysis.  Measurements  for  MPBX  2206  and  2207  apparently  indi¬ 
cate  the  development  of  a  tension  arch  and,  if  so,  this  would  not  have 
been  predicted  by  the  elastic  model  used  in  the  finite  element  analysis. 
Also,  the  direction  of  movement  probably  was  controlled  by  the  Joint  sys¬ 
tem  around  the  opening.  The  large  measured  displacement  shown  for  MPBX 
2217  probably  vas  caused  by  this  instrument  being  located  in  the  same 
large  wedge  of  moving  rock  as  MPBX  2215.  The  fact  that  the  direction  of 
the  motion  for  2217  is  opposite  to  that  for  2215  might  be  attributed  to 
the  possible  rotation  of  the  block  during  its  movement.  In  either  case 
a  finite  element  analysis  based  upon  an  elastic  model  would  not  be  ex¬ 
pected  to  predict  these  anomalies. 
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5.  CONCLUSIONS 


The  iDstrumentation 

1.  The  usefulness  of  the  MPBX  measurements  to  some  extent  depends  upon 
how  its  axis  is  oriented  with  respect  to  the  direction  of  Jointing, 
faults  and  shear  zones. 

2.  A  reasonable  assumption  is  that  a  considerable  release  of  elastic 
strain  occurs  immediately  upon  excavation  of  an  opening;  therefore,  the 
measuring  instmment  should  be  installed  as  rapidly  as  possible  after  the 
excavation  of  a  tunnel  section  is  completed. 

3.  Displacement  data  probably  can  be  lost  because  of  the  time  lag  between 
the  excavation  process  process  and  the  installation  of  the  Instrument; 
therefore,  it  would  be  desirable  to  install  strain-or  displacement¬ 
measuring  instruments  ahead  of  the  excavation.  To  minimize  construction 
disturbances  to  the  instrument,  the  most  desirable  installation  method 
would  be  to  place  the  instrument  in  a  hole  bored  from  the  ground  surface 
to  the  perimeter  of  the  proposed  excavation;  then  the  behavior  of  the 
rock  system  can  be  recorded  as  the  excavation  approaches  and  then  passes 
the  instrument. 

k.  It  is  possible  that  significant  but  unrecorded  strains  may  occur  in 
the  intervals  between  readings.  To  minimize  the  loss  of  such  data,  it 
would  be  desirable  to  attach  continuous  recorders  to  each  instrument.  If 
this  does  not  prove  feasible,  it  would  be  worthwhile  to  explore  the  devel¬ 
opment  of  a  "meucimum-miniraum"  strain  recorder  that  would  be  permanently  at¬ 
tached  to  each  cantilever  on  the  MPBX;  it  would  act  in  a  manner  somewhat 
analogous  to  a  maximum-minimum  thermometer,  i.e.,  it  would  be  designed  to 
indicate  the  maximum  and  the  mininaim  strains  that  occurred  in  the  time 
intervals  between  periodic  readings  of  the  cantilevers. 


Preceding  page  blank 
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Comparative  Analysis 

5.  The  results  from  our  study  indicate  that  a  finite  element  method  of 
analysis  that  is  based  upon  an  elastic  continuum  tends  to  predict  dis¬ 
placements  that  can  be  one  or  more  orders  of  magnitude  less  than  vhat 
actually  occurs,  and  occasionally  is  larger. 

6.  The  preceding  conclusion  appears  to  particularly  apply  when  the  rock 
system  contains  numerous  Joints  and/or  faults,  and  secondly  when  the 
rock  properties  used  as  input  to  the  analysis  are  derived  from  intact 
laboratory  specimens  rather  than  from  in  situ  field  tests. 

7.  When  the  rock  system  is  relatively  homogeneous  from  an  elastic  view¬ 
point  and  contains  few  Joints  and  no  faults,  the  finite  element  emalysis 
generally  but  not  always  can  predict  the  direction  of  the  displacements 
that  may  occur  during  excavation  although  there  still  may  be  discrepancies 
between  the  magnitudes  of  the  measured  and  the  predicted  displacements. 

8.  In  general,  since  our  study  indicates  that  the  accuracy  of  the  finite 
element  prediction  for  underground  openings  appears  to  be  highly  sensitive 
to  the  geological  environment,  finite  element  analyses  based  upon  elastic 
models  of  such  openings  should  be  used  with  considerable  caution.  The  pre¬ 
dictive  accuracy  of  the  finite  element  method  should  be  significantly  im¬ 
proved  if  (1),  a  method  can  be  found  to  quantify  all  significant  geological 
factors  so  they  can  be  used  as  direct  inputs  in  the  finite  element  model, 
and  (2)  data  are  avadlable  on  the  orientations  and  magnitudes  of  the  ambient 
stresses  around  the  opening. 
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Figure  11.  Straight  Creek  Tunnel  -  Sta.  78  +  17  -  Displacements 
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Figxire  15.  Straight  Creek  Tunnel  -  Sta.  114  +  53  -  Displacements 


^  tTNIHI  m4  0I»  MUCTS 
^  rAULT  f  SHIAH  ZOttt 
MOKIN  ZOHt 

^  CONTACT  MTWCIN  LlTHOUMHI 
w(l  MIASUNINI  STATION 
NO.  At  STATION 


T  SlOo 

1 

....  J 

_ A _ 1 

SCALC 

•  ^  ^  O* 


NOTI  : 

tvnntl  It  twrvtA  SaIvma  SfAtiM  Ho.  itT  tnS  SioNoo  No.  It*. 
AT  sitilon  No.  >17  tlio  dirootion  of  INo  tonoti  li  NSf*  OO'St*  W. 
AT  STOTtoo  No.  lit  TNo  Okooffoo  •«  TNo  hmiroi  SSS*Ot'tl*‘«. 


M  ^ 

Figure  16.  Pacheco  Tunnel  -  Geologic  Log 
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18.  Pacheco  Tunnel  -  Sta.  132  +  50  -  Displacements 
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Figure  19.  Pacheco  Tunnel  -  Sta.  151  +68.1  Displacements 
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Fig\ire  20.  Navajo  Tunnel  -  Geology 


Figure  22.  Navajo  Tunnel  -  sta.  7  +  46  -  Displacements 
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Figure  24.  Navajo 
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Figure  27.  Green  River  North  T\innel 


Figure  28.  Green  River  North  Tunnel  -  Sta.  103  +  51  -  Displacements 


Figxire  29.  Green  River  South  Tunnel  -  Sta.  105  +  05  -  Displacements 


Figxure  30.  Green  River  North  Tunnel  -  Sta.  105  +  54  -  Displacements 
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APPENDIX  A 


TIME-DISPLACEMENT  PLOTS 


Explanation; 

SC  -  Straight  Creek  Pilot  Bore 

GR  -  Green  River  Highway  Tunnels 

MP  -  Morrow  Point  Undergroxmd  Powerhouse  Chamber 

N  -  Navajo  Tunnel 

P  -  Pacheco  Tunnel 

Example : 

fs^  -|2001|  fllU  ♦  33 1 

^MPBX  Location  in  the  tunnel 

^  Number  of  the  MPBX 

*— Name  of  tunnel 

The  small  numbers  on  the  right  side  of  each  figure  refer 
to  the  suichor  represented  by  that  particular  plotted  line. 
The  reference  anchor  that  is  assumed  not  to  move  is  the 
horizontal  line  that  separates  the  (+)  and  the  (-) 
displacements . 
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Computer  Programs  for  Time-Displacement  Plots 

by 

Patricia  Nahas 


The  gOEil  of  the  computer  programming  aspect  of  this  research  on  strain 
distribution  vaa  the  production  of  plots  of  displacement  vs  time  around  each 
of  five  underground  openings.  (See  page  B-5  for  sample  printouts  of  the 
raw  data.)  Pvirdue  University's  CDC  6500  computer  and  CALCOMP  Plotting  System 
were  used  for  the  analyses.  The  CALCOMP  system  at  Purdue  consists  of  a  set 
of  subroutines  that  produce  a  file  on  a  disk  that  is  copied  onto  a  plot  tape. 
This  tape  is  then  run  on  a  Nfodel  563  CALCOMP  Digital  Increment  Plotter, 
producing  an  X-Y  plot  of  digital  data.  (The  user's  manual  for  the  CALCOMP 
Plotter  can  be  obtained  upon  request . ) 

The  Bureau  of  Reclamation  was  the  source  of  the  raw  data  for  the  Morrow 
Point,  Pacheco,  and  Navajo  Tunnels.  It  was  also  the  source  of  Program  111 
in  the  appendix.  This  program  produced  output  from  raw  data  that  included 
displacements  and  associated  times  for  each  anchor  for  each  station.  With 
slight  modifications,  this  program  was  run  on  the  CDC  computer  at  Purdue, 
and  in  addition  the  above-mentioned  portion  of  the  output  was  placed  on 
punched  cards.  (Program  #1  in  the  appendix  is  the  modified  version  of  the 
program  we  received  from  the  Bureau  of  Reclamation.)  The  data  on  the  result¬ 
ing  set  of  cards  were  exactly  the  points  to  be  plotted  from  the  strain 
distribution  curves;  Program  #2  in  the  appendix  is  the  program  used  to 
produce  the  plots. 

The  initial  part  of  the  routine  reads  in  the  identification  information 
and  stores  the  times  and  displacements  in  arrays.  The  next  section  of  the 
program  determines  the  maximum  and  minimum  displacements  and  this  information 
is  used  to  scale  the  Y  axis.  The  scale  for  the  X  or  time  axis  is  predetermined: 
it  can  be  uniform  for  all  the  data  or  it  can  be  set  to  reflect  the  time 
spread  for  the  particuleu:  data  in  any  one  set  of  plots.  In  any  event,  it 
is  necessary  to  know  the  time  span,  say  in  hours,  over  which  the  raw  readings 
were  collected.  The  last  part  of  the  program  consists  mainly  of  calls  to 
the  subroutines  which  produce  the  plots.  These  calls  provide  information 
regarding  the  length  and  labelling  of  the  axes  and  the  points  to  be  plotted. 
Finally,  the  program  tests  the  counter  variable  KOTD  against  the  constant 
NOT  which  is  set  at  the  start  of  the  program.  (Note  that  NOT  could  be  read 
in  for  greater  flexibility.)  NOT  indicates  the  number  of  data  sets  -  and 
hence  the  number  of  plots  -  to  be  processed  on  a  peurticular  run  of  the  program. 
If  NOTD  does  not  equal  NOT,  the  pen  must  be  repositioned  before  another  plot 
is  begun.  Otherwise,  the  program  is  finished. 

In  the  case  of  the  data  on  the  Green  River  and  Straight  Creek  Tunnels, 
it  was  necessary  to  write  the  routine  to  compute  the  times  and  displacements 
from  the  raw  data.  Programs  113  and  0k  in  the  appendix  are  the  programs  used 
to  convert  the  data  from  Straight  Creek  and  Green  River  respectively.  These 
progi'ams  also  produce  the  plots  for  the  data  being  processed.  They  are 
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essentially  alike,  varying  only  where  necessary  to  accommodate  the  slightly 
different  formats  of  the  data  and  the  different  labellings  of  the  plots. 

The  initied  part  of  the  programs  reads  in  identification  information, 
the  calibration  factor  for  each  anchor,  and  the  raw  data.  Taking  the  time 
of  the  first  reading  as  the  zero  reference  time,  the  program  then  determines 
the  time  in  hours  of  each  successive  reading.  The  following  section  of  the 
program  is  devoted  to  the  computation  of  the  displacements,  with  those  of 
the  reference  anchor  necessarily  first.  It  also  includes  procedures  for 
eliminating  bad  readings  from  the  raw  data*,  i.e.,  readings  of  00000.  or  99999 • 

In  general,  the  policy  with  respect  to  bad  readings  was  as  follows. 

If  the  bad  reading  occurred  for  the  reference  anchor,  it  was  eliminated  and 
along  with  it  all  corresponding  readings  for  the  remaining  anchors,  even  if 
they  happened  to  be  good.  Hence,  referring  in  the  appendix  to  the  raw  data 
for  Straight  Creek  Tunnel  (page  B-  5  ) ,  lOl*  +  83,  #200U,  note  that  three 
readings  were  eliminated  for  anchor  (i'S,  and  correspondingly  for  the  other 
seven  anchors,  even  though  the  readings  for  anchors  #2,  #3,  and  ffh  were  not 
all  bad.  In  the  case  where  a  bad  reading  occurred  for  an  anchor  and  the 
corresponding  reading  for  the  reference  anchor  was  good,  the  reading  was 
eliminated  only  for  the  anchor  in  question.  Thus  with  reference  to  the  same 
raw  data  cited  above,  one  counts  2h  good  readings  for  the  reference  anchor, 
but  only  17  for  anchor  1/1 ,  23  for  anchors  02  and  #3,  19  for  anchor 
15  for  anchor  #5,  20  for  anchor  #6,  and  l6  for  anchor  #7.  This  procedure 
for  eliminating  bad  data  was  selected  from  among  several  alternatives. 

Although  it  complicates  the  programming  -  since  one  must  keep  track  of  the 
differing  number  of  points  to  be  plotted  for  each  anchor  while  medntaining 
the  correct  time  value  across  all  anchors  -  it  enabled  us  to  use  as  much 
of  the  available  data  as  possible. 

An  important  note  of  information  is  that  the  readings  for  the  reference 
anchor  are  tested  first  and  if  they  are  bad,  they  are  eliminated  along  with 
the  corresponding  readings  for  the  other  anchors.  The  time  array  is  edso 
altered  accordingly  and  the  counter  NOR  -  the  variable  indicating  the  total 
number  of  readings  for  the  data  set  in  question  -  is  decreased.  If  the 
reading  is  good,  the  program  computes  the  displacement  for  the  reference 
anchor  as  follows : 

D(K,NA)  =  [R(K,NA)  -  R(1,NA)]  *  CAL(NA) 

where 

D(K,NA)  =  displacement  for  the  reference  anchor 
R(K,NA)  =  reading  for  the  reference  anchor 
CAL(NA)  =  calibration  factor  for  the  reference  anchor 

Once  the  displacements  for  the  reference  anchor  have  been  determined,  the 
program  proceeds  to  successively  compute  those  of  the  remaining  anchors.  It 
must  take  into  account  the  procedure  for  eliminating  bad  readings.  This 
involves  maintaining  a  count  of  the  number  of  good  readings  for  the  anchor 
in  question  and  use  of  the  correct  readings  for  the  reference  anchor  in  the 
computations  of  the  displacement.  The  general  formula  for  the  displacement 

D(K,J)»  [R(K,NA)  -  R(1,NA)]  •  CAL(NA)  -  [R(K,J)  -  R(1,J)]  *  CAL(j) 
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where 


D(K,J)  ■  displacement  for  the  anchor 

R(K,J)  =  reading  for  the  anchor 

CAL(j)  =  calibration  factor  for  the  anchor 
til 

R(K,NA)  =  K  reading  for  the  reference  anchor 
CAL(NA)  =  calibration  factor  for  the  reference  anchor 
emd  D(1,J)  *  0  for  every  J 

However,  aince  the  array  of  raw  readings  for  a  particular  anchor  is  corrected 
before  the  displacements  for  that  anchor  are  computed,  the  n^^  reading  for 
anchor  J  may,  for  example,  correspond  timewise  to  the  reading  for  the 
reference  anchor.  In  this  case,  the  general  formula  for  computing  the  dis¬ 
placement  must  be  altered  and  appears  as  in  statement  with  the  correction 
variable  KOUNT.  Reference  once  more  to  the  case  of  Straight  Creek  Tunnel, 
lOU  +  83,  #200J*  may  help  clarify  the  situation.  We  have 

D(l,8)  =  0 

D(2,8)  =  (  10050  -  13535  ) 

D(3,8)  =  (  10000  -  13535  ) 

D(l*,8)  »  (  1U570  -  13535  ) 

1^5,8)  =  (  10060  -  13535  ) 

D(6,8)  «  (  1U735  -  13535  ) 

and 

D(l,l)  =  0 

D(2,l)  =  D(2,8)  -  (  11920  -  12975  )  •  (  112  ) 

D(3,l)  =  D(6,8)  -  (  1U2OO  -  12975  )  •  (  112  ) 

=  (  IU735  -  13535  )  •  (  106  )  -  (  IU2OO  -  12975  )  .  (  112  ) 

where  D(3,1)  corresponds  to 

D(K,J)  =  [R(K+K0UN1',NA)  -  R(1,NA)]  •  CAL(KA)  -  [R(K,J)  -  R(l,J)]  •  CAL(j) 
with 


•  (  106  ) 

•  (  106  ) 

•  (  106  ) 

•  {  106  ) 
•  (  106  ) 


K  *  3  NA  =8 

J  =  1  KOUHT  =  3 

The  next  section  of  the  program  determines  the  maximum  Euid  miuimijn 
displacements  and  thus  determines  a  seeding  factor  for  the  plots.  The  times 
and  displacements  to  be  plotted  for  each  anchor  are  then  printed  out.  These 
values  do  not  appear  explicitly  on  the  plots  and  a  listing  of  them  enables 
one  to  know  exactly  the  X-Y  coordinates  of  a  particular  point  on  the  plot. 
This  listing  can  also  be  used  to  check  that  the  computational  portion  of  the 
routine  is  correct. 

The  latter  part  of  the  program  is  devoted  to  the  plotting.  Calls  to 
various  subroutines  provide  the  information  necessary  to  draw  and  label  the 
axes  and  plots.  Since  only  one  line  can  be  plotted  at  a  time  but  plots  for 
all  anchors  on  one  set  of  axes  are  desired,  the  points  to  be  plotted,  stored 
originally  in  two-dimensioned  arrays,  are  tremsferred  to  one-dimensional 
arrays  inside  the  DO-LOOP  that  actually  plots  the  line.  This  fact  plus  the 
fact  that  the  last  two  elements  of  the  one-dimensional  arrays  must  provide 
information  on  the  number  of  points  to  be  plotted  and  the  scaling  factor 
indicate  why  such  care  was  taken  when  the  displacements  for  each  anchor  were 
determined. 
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Lastly,  the  program  checks  to  verify  that  all  plots  for  a  particular 
run  have  been  processed  and  either  ends  or  returns  to  read  the  next  set  of 
raw  data. 

The  above  hopefully  will  prove  helpful  to  those  who  wish  to  repeat  these 
procedures  as  part  of  their  own  research.  It  is  not  possible  to  go  into 
complete  detail  in  describing  the  data  and  programs  used.  Some  knowledge  of 
computer  techniques  is  assumed,  and  it  is  believed  that  by  reference  to  the 
appendix,  the  reader  can  reproduce  the  method  used  and  understand  the 
rationale  behind  it. 
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Example  of  haw  Strain  Data  irom  Straight  Creek  Tunnel 


00000? 

00000? 

00000? 


OOOOO? 

00000? 

000014 

000014 

0000?? 

0000?? 

00003? 

00003? 

000033 

000034 

00003S 

00Q03A 

0000<^1 

000061 

O0QO64 

000066 

000066 

000070 

000071 

00007? 

000073 

000076 

OOOIOS 

000113 

OOOlU 

OOOllS 

00011S 

000117 

0001?? 

00013? 

000134 

000136 

000140 

000141 

00014? 

000146 

00O161 

00016? 

000166 

000166 

O0016O 

000164 

000164 

000166 

000174 

000176 

000176 

000177 


COMPUTER  PROGRAM  NO.  1 

pqOfiBAM  EXPEH(IwPUT,OUTPUT#TAPE5«lNPUTtT4PE6»0UTPUT«PUNCH.T4PE7  ■ 

IPIJNCH) 

C  MPPX  DATA  REOUCTION  PPOfiRAM  NUMERTr  AND  CAROS  FOR  X  V  PLOT  W  ORTEL 

niMENSION 

OTMFNSION  UTSPL ffl«300) »FT(300) 

DTMFNSION  A0FPTm(8)  ,OCAnT  (b*300 )  iRrANT  (8,300)  *DAV(300)  fAMONTHOOO) 
1,YFAR(300)  iHOUROOO)  fFEAT(40)  iDELTa  (R*30h)  *nAYM0(12)  #SUM(8#300)  *TI 
?mF(300) fDELTAb (fl,300) #OlFF(8»300) «bElTA (8,300) »BELTA8 (8,300) (Ani^S 
3T(B,100) 

COMMON  TIMEnOOo)  tRnG(b,lO)  *Cl  (8,16)  ,C2(8»10)  iN,NROQtHCANT 
S  PFAn(S»10)  (FEAT (I) *1*1*40) tSTAtPRvM 
10  fORmAT(40A1,?AR) 
pFA0(5»B)  MCANT 
R  format (14) 

PFAn{5*?0)  (AOEpTH(T) •I*!*^) ‘C 
?0  format (flFb.O^FS, 4) 

CALL  CALRDG 
Nal 

rltlst  *  0*0 
nLTl  ST  *  0,0 

?6  RFAn(St30)  (UCAnT(I*N) tlaltB) •AMONTH(N) *nAV(N) »VEAR(N) fHOUR(N) *  (RC 
IANT(I,N) ,Isl«8) 

30  FORMAT(BFt>.0*lX,3F2.0*FS*2»20X  /BF6.0*3?y) 

TF(nAY(N)^AMONTw(N) )  40,40*35 
36  NsN*l 

00  TO  ?5 
40  NsN*1 

no  60  I«l«8 

KsO 

no  60  J*1»N 

TF(0CANT(I*J)  )  89*50*53 

63  nTFF(I*J)«(UCANT(I*J)-HCANT(I*J)  )/?,0 
T(f(y)  55*55,64 

64  04Y0IF*0IFF(I*J)-DIFF(I*K) 

PELTA(I*J)»UAYDTF 

Tl*T 

JJ*J 

CALL  CALNO(DCANt,II,JJ*Y1) 

CALL  CALNO(0CANt,II,K*Y2) 

Y3«AflS(Yl-Y2) 

XCANTsABSIOCANT (7* J)-0CANT(I*K) ) 

TF(XCANT*EQ,0.)  60  TO  72 
FACT*(Y3/XCANf ) 

60  TO  ^5 
72  FACT  ■  0* 

75  FNUMH-F act*  (1*0* {C*40EPTH ( 1 )  )  ) 
nFLTAd  ,J)*L)AYOTF*FMUMb 
GO  TO  56 

69  0FLTA(I,J)«DLTL6T 
RELTA(I*J)«8LTL6T 
GO  TO  60 

65  nFLTA(I,J)*0.0 
RELTA(I,J)*0.V 
GO  TO  57 

66  nFLTA(l,J)«OELTA(I, J)4UELTA(I*K) 

RFLTA{I*J)*BELTa (I*J)*BELTA(I*K) 
nLTLST«l)ELTA(I,j) 

PLTLST«BELTA(I, j) 


00200 

:oo?o? 

00207 

00210 

j00277 

00221 

[0022S 

’002?7 

00240 

00241 

00241 

0024S 

00247 

00247 

|00247 

|o0247 

e002'SI 

002SI 

002S4 

002*»6 

002^0 

002AA 

002^*7 

00271 

0027^ 

00300 

0030? 

00304 

(00311 

{00317 

00320 

0032? 

[00320 

S003?6 

100330 

<00349 


0034? 

00341 

00346 

00310 

100312 

<00377 

*00377 

100401 

•  00403 

•  00406 
(00424 


•00494 


17 

60 

601 

104 

101 

102 

106 

103 

101 


62 

61 


70 


KSvi 

CONTINUE 

reap (5*501  )  i*L. (An JUST (I*L) *I«l*6* 
rORMAT(2I3*flF«.6) 

TF(K)  61t6ltl02 
J*L 

no  103  I«l*fl  ,  , 

nELT A ( I • J) *OELT  4  ( 1 1 J) ♦An JUST ( I *L ) 

CONTINUE 

IF(N-J)  504t104«105 
J*J*  I 
C,0  TO  506 

nATA^^UAYMO ( T ) «T*i,]2)/0.0*3l#0*2«.0*31,ft*30.0*31#0*30«0*31,0»3l. 

I0*30«0»31«0»30,o/ 
no  10  I*ltN 
A*0«0  ^ 

TF(VEAR<I)-YEAPfl) »  62»65f62 

4*361.0 
L*AM0NTM(I) 
no  70  J*1*L 
A*A^0AYM0(J) 
continue 

TIMF0(1)«TIME  (n-TlME(l) 

60  continue 
no  61  1*1*8 
no  65  J«l*N 

1UM(I*J)*0CANT  (t»J)*RCANT {l*J) 

0FLTA8  ( I  ♦  J)  *UELt A  ( I ,  J) 

PELTASd  *J)*NtLTA(I»J)-6ELlA(MCANT,  J) 

CONTINUE 
MsO 

uoTTF(6a1001  ( FFAT ( T  1  *  I *1  * AO)  * ST A  * PRXM 

rOR«ATtlH?.*0*.»OHMPB«  e«TtN»0METEO  //mxJI 

1 i A* 1 / // . AlX. TMSTATlONt Afl.6Xf 8MWPBX  MO,  ♦  A6//»53X*8HDATF  //.53X  l 

i  ,  cant,  VR  E  C.RTUVR  3  C.RTLVR  * 

3CANTLVR  5  CAnTLVP  6  CANTLVM  7  CANTi  VR  8  DATE  TyME  //) 

HaH*6 

10AV*0AY (K) 
month  ■  AM0NTH(«) 

IS“E‘tI|!o!''lOrANT(I.K..IM.8I.MnMTH.ID.».IYEAR.HOOR(K..(RfANTlI. 

iK) *  1*1 *8) * (SUN  (  y  ,K) f I*l *8)  T9  Fa  9  /  16HRE 

?00  format  (16H0I9ECT  REAOlNOS  10*  '  *■***  2,yH/,I2*lH/*  I2,F6.2  /.16HRE 
IVFRSE  RFAnlNfi^*8F10.1^»l5HSUM  (^F  RfAOING6»6F10»1//) 

KsK*1 

IF(N-K)  250t2l0,210 

iJcE  DiIpl°F»OM  N0,y2*30HOlSPL  FROm  FACE  OISPL  FROM  mO.I?/, 56**60 
♦h(IMD  UNITS)  (TNO  UNITS)  (INCHES)  (InCheD/) 

Ma3 


61 


60 

100 


100 


B7 


0004?S 

0004?7 

000431 

000433 

0004^0 

0004^1 

0004S? 

0004S4 

0004S^ 

0004S7 

000473 

000473 

00047^ 

0005?? 

0005?? 

000527 

00053? 

00053? 

000534 

000536 

000544 

000546 

000551 

00055? 

00056? 

00056? 

000565 

000603 

000605 

000606 

000607 

000611 

000625 

000626 

000633 

000641 

000645 

000647 

000651 

000663 

000663 

000666 

000666 

PROGRAM 

070306 


no  475  J«1*N  i 

TF(M-J)  47lf47l,460  ’ 

471  M«M*3 

WRlTE(6t300)  (FFAT(T)*lal*40) •STA,PRXHfMCANT«MCANT 
460  rONTINUE 
inAY*0AY(J) 

MONTHbAHONTH ( J) 

IYEARbYEARU> 

ITIME«TIME0( J) 

WRITE(6«4<!5)  MOMTHtTOAY,IYEARiITlMr 
425  F0RMAT(2Xi4HOATF*l2.1H/tI2»lH/#I2/,2¥,l6HELAPSE0  TIME <MR) , 15/) 

no  475  i«i»a  '  ^ 

445  WRITE<6t450)  I  * aOEPTH(I ) #BELTA ( I » Ji ,RELTA8 ( I » J)  •DELTA ( I . J) .OELTAG ( 
1I«J) 

450  FORMAT<29X.I2.1?X#F5,l#10X.F7.0#10x,r7.n,8X,F10,6t7X,FlO,6) 

475  CONTINUE 

WRITE(6«500) 

500  format (IMI) 

00  591  Ib1«N 
FT(I)  ■  TIMEO(I) 

00  591  Kal«a 

OlSPL(K,I)  a  OElTA8(K»I) 

591  continue 

00  100«  Ial*N 

PUNCH  4.FT(1) . (nlSPL(KtI)tKal,8) 

4  F0RMAT(1X«F6.0,RF9,6) 

1000  CONTINUE 

00  600  I«1*N 

STRAIN ( 1 . 1 ) a (-1 ,0) aoELTA  <MCANT» I ) -OELTAR  J 1*1) 
STRAlN(l,I)aSTRAlN(l,I)/(ADEPTH(I)#l?,0) 

600  CONTINUE 

00  700  Ia2,MCANT 
00  700  JaltN 
I2aT-l 

STRAIN ( I  •  J )  a ( 0ELTA8 1 12* J) -OELTAB ( I , J) ) /  (  (  ADEPTH ( I ) -ADEPTM ( I ? )  )  • 

11?. 0) 

700  CONTINUE 

WRITE(6*S00) 
no  aoo  lal.N 
ITlMEaTIMEO(I) 

WRITE(6«900)  I.TTIMFt(STRAlN(J»I)  , jalfB) 

900  FORMAT(l5Xtl5|?K,i5,8FlO*6) 
aoo  CONTINUE 
00  TO  5 
FNO 

I  : 

LENGTH  INCLUDING  I/O  BUFFERS 


UNIISFD  compiler  SPACE 
053700 


0000^1 
Qooom 
000001 
000003 
I  000010 
000010 
'  00001? 

0000?4 
I  0000?4 
'  000031 
00003? 
i 000040 
I  000041 
f 000046 

t  0000^0 

00005^ 

0000S6 


5 


8 


SUBROUTINE  CALRnO 

OTMFNSION  OEPUTn(8.10>  .  ^  . 

COMMON  TIMEO(300) fRnO(8«10)  *C1 (8*lft) tCECRt 

no  10  i«i#0 
read (5*5)  NR06 
format (!♦) 
no  10  J«l*NR06 

RFAn(5*8)  OEFLTn(!»J) *BoO(I*J) 
format (2F8.0*64k) 


10) vNcNROOtMCANT 


10  continue 
no  ?o  i»i*8 
no  ?0  J»2#NR06 

c“  n  1k)  -  (OEFLTN  ( I .  J)  -OEFLTN  ( I  *K) )  /  » R06  ( 1 ,  J)  -RDO  ( I  .K)  ) 

C?(1*K)«0EFLTN(T»J)-C1  (I*F)*R06(f  *.)) 


20  CONTINUE 
return 


FNO 


iSURPRORRAM  length 


i  000??3 

UNUSFn  compiler  space 

lOSTOOO 


000006 

000006 

000006 

000017 

0000?! 

0000?? 

0000?3 

000035 

000036 


SUBROUTINE  CALNO  (X.I.J.V) 

CnMMON*TIMEO(3OO)*RnG(0*lO)  tCKB*!")  *02(8*10)  *N*NROG*mCANT 

no  5  k«i*nrog 

IF(R00(I*K)-X(I*J))  6*6*5 
6  MAXbK 
S  CONTINUE 

V*C? ( I *MAX ) ♦Cl  ( T  *MAX ) *X ( I  * J) 

return 

FND 


subprogram  length 
000062 


UNUSFn  COMPILER  SPACE 
057000 


COMPUTER  PROGRAM  NO.  2 


UUb  ub  i! 

Ub  2 
li  0  u  j  C  ^ 
jbO  0  b 
C  Uu  U  t  3 
w  wO  J  1  9 
bOOuXS 

J  Ju  li  2  7 
b0bG4X 
w  bOu  41 
ObbUH3 
OUG  b  3  7 
u  ju  b  37 
bUtj  Gbu 
GGUGb2 

bOUbbi 
b  U  U  b  /  U 

U  Ub  U  7  b 
u  bO  b  7  3 

UUb lb  U 

C  bb  1  b  2 
QUb lb  3 
juOiGh 
duullG 
auuiii 

b  0 J 11 4 

0aul21 
aMul24 
uCibi3  J 
JUu1j92 
UOG133 
UUbl43 
0001^4 
u01il46 
uGb  147 
uiij  151 
aObl53 
ttuulbi 
biiu  lb5 
buwl7l 
ugbl75 
GObZCl 
00u2C5 

ki  00211 

0u02l7 

00u2o2 

bTffTsT 
000254 
uOu  257 
'000261 
u  bb27b 
a03272 
~JOuii/6 


PkuGRftM  RlCl( INPUT, OUT^JT, PLOT, T A >£5*1 NPOT,TAPE6»OUrPUTtTl»£isPL01 
1) 

OiliCN^iUN  FtLAr(HJ),TX(l5J>,R(8,l5j),0SY(13bi 
CALL  PlOTj 
NOT  *  <. 

NOTD  *  1 

Ibbw  RlAU(9,a)  NORfNOA 

1  FORMAT (214) 

RtA0(9,2)  MO,IOA,1YR 

2  FuRHAT (A3,Ij, 15) 

RLA0(9,j)  (F£AT(I)  ,I*l,4ti)  ,STA,P3<H 

3  F0RHAr(4jAl,2AOJ 
00  4  <  =  l,NUr< 

4  »<lA0(3,9)  T  X  (<)  ,  (R(J,K)  ,  J*1,N0A) 

5  FORMAT (AX,Fb. b,8F9.b) 

YMAX  =  R(l,l) 

JO  lb  b*l,NOA 
JO  xj  K=2,NOR 
IF{R(J,K) .GT. YMAX)  GO  TO  9 

ijO  T U  1 . 

9  YMAX  -  R(J,K) 
lb  JCNTlNUc. 

YMIN  =  K(x,x) 

30  li  J*1,NUA 
JO  13  K*2, NOR 

IF(R(J,<) .lT. YMIN)  GO  TO  1h 
GO  TO  l9 

14  VMIN  *  R(J,K) 

15  CONTiNUc 

YSCALE  *  (YMAX-YMIN)/9.b 

CALL  SYMBOL (4.  ,*.4,. 2, FEAT, 4.0,40 ) 

NNsNOR+l 
_NNN*NN*1 
30  2j  J*1,N0A 
R(J,NN)  *  b. 

2b  R(J,NNN)  sYSGAlE 
'TXTTTN)  i'oV 
TX(NNN)  *  133b. 4 

call  3YMBOL(1.0,10.3,.l4,8HSTATION-,b. ,0) 

CALL  SYMBOL (2.3,l0.3,.lb ,STA , 0 . o , 3) 

CALL  SYMBOL (3.8,10.3,.lb,9HMBPX  NJ.>,1.,9) 
call  SYM8OL(5.2,lb.3,.lO,P8XM,o.ii,0) 

call  SYMBOL (1.0,  9.9,. lb,29H0ATE  OF  FIRST  ACTIVE  READING-, 0 ., 29) 
CALL  SYMB0L(5.3,  9.9,.lo ,NO,u. ,3) 
call  NUMBER(3.a,  9.9,. iU , ZOA ,o .  , 2NI 3) 

CIILL  ■NLIHSERT6.3,  ^.9,.ib,I/R,0.,6^1M,,15) 

CALL  PLUT(o  .  ,  (  (-YMlN/YSCAL£)^.4b)  ,°3) 

CALL  AXlSib  .  ,  (YMIN/YSCAL£),22m0ISPLACENENT  IN  INCME5,4'22,9.  o,  9U.U, 
lYMlN^YSCALt,-!) 

CALL  AX1S(0. , (YMIN/YSCALE)  ,11HTIME  IN  HRS,-11, 15., 0. 0 ,TX (NN) , TX (NN 
lN),-u) 

CALX'  Plot  ( 0 .  . 

CALL  PlOT(15. ,0.*42) 

00  50  J*l,NOA 

00  40  K»1,'NNN  '  "  .  - 

40  OSY(<)  «  R(J,K> 

CALL  LlNE(TX,OSV,NOK,lfM,Ol 

eiLL  ~Wff(B]ER'<  1 5 .  •  TVn  ^TTSIT  ( NMID  FT*  1 0  •  J,  S  •  B  ,  ZHZ2 1  - 


OOj  3o  7 

7k 

continue 

u003lu 

CALL  ®LOT(l?.| ( (TMIN/Y SCALE) -.46)  ,-3) 

uuUilb 

XF(NOrO.cU.NUT)  GO  TO  ISu 

UuO'317 

NOTO  s  NOTO  *  1 

uuO  32  u 

(.0  T  0  lb  b  u 

uOu321 

13u 

CALL  S YH30L  ( J  .  I  k  .  f  . 1 y 7 9  y9u*yi) 

001326 

CALL  PlOT(u tJ»993) 

0Uu33u 

STOP 

0uu332 

END 

- -  “ 

- - 

— 

PROGRAH 

006642 

LENGTH 

INCLUDING  I/O  BUFFERS 

-  — —  -  — 

-  -  ■  ■-  •  -  • 

UNUSED 

OlSbCu 

COHPILER  SPACE 

TTr2WI3.  31711/7 1. PURDUE  H4  3c  7u/i2/j>. 


'01.24.^3.HP2&6V  J66T4|JU0OVCS6OOCu7fPiTT2'Ow  *P2 
ul*24«23 

01.24. 23. HAP(ON» 

TiTir.liiTFOHT^T -  - -  '  " 

U1.24.25.CTIHE  JuO.oui  SEC.  FUM  NOO  LEVEL  buE 
al.24.25.L60. 

01. 24. 26 .CX  2.U64  SLC.y  P%  3.2l?  fEC. 

Ul.24.26.NL  1620ti  NOROS 

Ul.24. 27. PLOTTING  STARTED 
ai.24.3l.PLOttlNG  S'JCCESSFUL' 

01.24. 31. STOP 

01.24.31.COPVPLT(PERLOFFfl6073fllX  6u, 2, BLACK) 
01.24. Sl.Cx  $.406  SEC.f  ~TX  "4.131  STC. 

01.24.31.nl  4uuO  NOROS 

02.47.09.  NT51  ASSIGNED  •  OALCOHP 
02.-47.5l.C0PVPLT  Cd^PieTEO— EXPECT  PLOT  OUTPUT 
u2. 47 . 61 .PLOT  0.3b7  HR. 

02.40.00. CP  0.792  SEC..  PP  42.930  SEC. 

flZ.U.uO.L'XKieS - 2^^^ -  - 

02.4O.u0.CH  j.2b4  NND-SEC.f  FL  4uuU  WORDS 


BU 


COMPUTER  PROGRAM  NO.  3 


PARTIAL 

SN 

NA 

NOR 

NOA 

CAL(J) 

DAY(K) 

HR(K) 

TIMIN(K) 

R(K,J) 

D(K,J) 

T(K,J) 

NOGR(J) 

NOT 

NOTD 


INDEX  TO  VARIABLE  NAMES  THAT  APPEAR  IN  PROGRAMS  ftZ  AND  #U 

Station  number  for  data  being  processed 

Reference  anchor  number  for  data  being  processed 

Number  of  readings 

Number  of  anchors 

Calibration  factor  for  anchor 

Day  of  reading 

Hovir  of  reading 

Minutes  of  reading 

reading  for  anchor 
K  displacement  for  J  anchor 
time  in  hours  for  anchor 
Number  of  good  readings  for  anchor 
Number  of  sets  of  data  to  be  processed 
Counter  which  keeps  track  of  nvunber  of  sets  of 
data  that  have  been  processed 


B12 


PflORRAH  JUOO  ( INPUT  t  OUTPUT •  HLOT*  TAP25b  I NPUT .  T*PE6«0UTPi )T «  TAPF  1  bPLOT 

1) 

DIMENSION  N06R(P)  tCAL  <8)  •OAY(UO)  *HR(1«0)  tTlMlNdAO)  tR(  l^OtH)  tT  ( 14 

lOOOO? 

iPtR)  tUlUOfB)  «06V(140)  •tX(14a) 

i 

lOOOO? 

CALI.  PLOTS 

% 

lOOAOl 

NOT  «  5 

1)00004 

NOTD  «  1 

f 

* 

[lOOOOS 

1000 

PEAD(b»l)  NA 

1 

>00013 

1 

format (12) 

>00013 

PFAD(5«2)  NORfNOA 

00023 

2 

F0RMAT(/I3»I2) 

i 

00023 

RF*0(St3)  SN 

00031 

3 

FORmAT(29X»F5,0) 

00031 

REAr)(6t4)  (CAL(.J)  tJalfNOA) 

* 

00040 

4 

F0RMAT(/8XfflF6.?) 

00040 

no  6  k«i»nor 

00047 

6 

PEAn(5t6)  OAT(K)  tHR(K) vT IMIN (K) t (» (Kt J) •J»l»NOA) 

00067 

6 

FORMAT(3F4*0<8FA.?) 

00067 

T2ERO  ■  24.*0AY(1)  *  MR(1)  ♦  TIMJN(n/60. 

00074 

DO  7  JaltNOA 

00103 

7 

T(I«J)  >  0*0 

00106 

no  g  kb2»noo 

00106 

no  6  JaltNOA 

00127 

6 

T(K»J)  ■  2«.«OAV(K)  ♦  HR(K)  ♦  TIMTN(K) /60.-T?FRO 

00131 

g 

CONTINUE 

00133 

LL  ■  NA-1 

00136 

no  10  jBltNA 

00146 

10 

ndfj)  c  0.0 

00147 

no  16  KaZtNOR 

00150 

11 

IF((R(K#NA),E0.gg99g.),0R.(H(K*NA),E0. 00000.))  60  TO  12 

00162 

n(KtNA)  M  (R(KfNA)-R(l  tNA)  )*;:AL  (NA) 

0O177. 

60  TO  15 

00172 

1? 

TF(K.EO.NOR)  go  to  14 

00174 

nor  ■  NOR-1 

00175 

no  13  Ibk.nor 

0020? 

R(I,NA)  ■  R(I«1,NA) 

00203 

T(I»NA)  «  T(I*l,NA) 

00205 

no  13  LalfLL 

00216 

R(I.L)  ■  R(I*lfL) 

1 

00216 

13 

T(ItL)  ■  T(I*1.L) 

00223 

60  TO  11 

00223 

14 

nor  b  NOR-1 

00226 

15 

NOOR(NA)  B  nor 

] 

00227 

16 

CONTINUE 

j 

100232 

no  17  Kb1,LL 

1 

00237 

17 

NOGR(K)  B  NOGR(NA) 

j 

00241 

KNOA  b  NOA-1 

1 

00243 

DO  60  Jb It KNOA 

00245 

MOUNT  ■  0 

00245 

NRR  8  NOOR(J) 

00247 

no  50  K82,NBB 

00251 

KKEEP  ■  K 

00252 

IF((R(KtJ), EO, 99999,). OR. (M(KtJ), EO. OOOPO, )  )  60  TO  20 

00263 

IF(KOUNT.OT,0)  GO  TO  40 

00266 

n(KtJ)  B-(R(KtJ)-R(l,Jh«CALcJ)*n{K.NA) 

00300 

60  TO  50 

00300 

20 

IF(KKEEP,EQ.MOOR(J) )  GO  TO  4^ 

00303 

25 

MOUNT  B  KOUNT*! 

B13 

1 

-j 

i 

'i 

j 

1  1  iidii 

000306 

000310 

000311 

00O32? 

000323 

00032S 

000336 

000360 

000361 

000363 

000364 

000371 

000373 

000374 

000376 

000377 

00040S 

000406 

000411 

000416 

000420 

000421 

000423 

000424 

000431 

000431 

000436 

000442 

000446 

000447 

000456 

000456 

000460 

000462 

000600 

000500 

000506 

000521 

000523 

000523 

000624 

000526 

000526 

000533 

000537 

000543 

000546 

000547 

000554 

000567 

000606 

000611 

000614 

000617 

000620 


30  N06R(J)  «  NOGR(J)-l  i 

NRR  ■  NOGR(j; 
on  35  KK«K6EEP»NBR 

0(KKf  J)  *  R(Ki^41  , J)  I 

35  T»KK*J)  »  r(KK*1»J)  I 

IF((R(KKEEP»J),FO,RR999,),OH,(R(KKEEP*J),EO.OOOOO.))  60  TO  20  | 

40  n(K«J)  ■-(r<(KfJ)-R(l*J))*C6L(J)*(R{K*K0l»igT*NA)-R(l,NA))*CAL(NA)  i 

60  TO  50  I 

45  N06R(J)  ■  N06R(J)-l  , 

MRR  m  NOGR(J) 

50  CONTINUE  | 

YMAX  ■  Odtl) 
no  200  JalfNOA 
NBR  >  NU6R(J) 

no  200  I«2«NflP  I 

IF(n(I,J),GT.YM4X)  GO  TO  IV9  ! 

GO  TO  200  ! 

199  YMAX  m  0(IiJ) 

200  CONTINUE 
YMIN  «  0(1*1) 
no  202  J«l.NOA 
NRR  ■  NOGR(J) 
no  202  I*2*NHR 

IF(D(I*J).LT.YMTN)  60  TO  201 
60  TO  202 

201  YMIN  m  U(IiJ) 

202  CONTINUE 

YSC4LE  ■  (YMAX-YMIN)/9, 
no  700  J«i»N0A 
WRTTE(6»750)  J*N06R(J) 

750  F0RMAT(1MI»4HJ  •  ♦I2«10X#23Mm0,  OF  GOOD  READINGS  •  ,13)  i 

NRR  ■  NOGR(J)  I 

no  700  KslfNBR 

WRITE(6*751)  R(K*J)  vTIKtU)  *D(K*J) 

751  FORMATdM  flOHRpAOlNG  «  ,F10 , 2, 1  OX , 7HT IME  ■  *F10.2*10X*l5HniSP|  AC 
IWENT  «  *F20,2) 

700  CONTINUE  , 

no  75  J»l*N0A 
NN«N0GR(J)*1 
NNN«NN«1 
n(NN,J)  ■  0. 
n(NNN*  J)  ■  yscai  e 
T(NN,J)  m  C, 

75  T(NNN,-J)  •  170.0 

CALL  SrMBOL(3,,l0.3,.15,9MNPBX  Nn.-,0.*9) 

CALL  NUMBER{4,3,l0.3f.l5fSN,0.,4HF5.0) 

KN  ■  NOGR(l)  «1 
KNN  ■  N0GR(1)  «2 

CALL  PLOT(0,, ( (-YMIN/TSCALE) *,46) .-3) 

CALL  AXIS(0.*  (YMIN/ySCALE),24M0ISPLACEMENT  TN  HTCRO-!N**24,9,0.90 
lO.YMINfYStALef-l) 

CALL  AXIS(0.*(YNIN/YSCALE),11HTIME  IN  HRS,-1  i  ,  15,0#0.0*T  (KN,  1 )  *T  d 
lNN,I),-0) 

CALL  PLOT(0.*0.,*3) 

CALL  PLOT(15.,0.,«2) 

00  100  Jb1,nOA 
MNbNOGR( J) 

NOGRJbNOOR(J) *2 


Bl** 


000622 

DO  80  KsltNOQRJ 

000632 

DSY(K)  *  D(K«J) 

000633 

80 

TX<X)  ■  T(X*J) 

000635 

CALL  LINE(TX«OSY«MN*ltOfO) 

000641 

CALL  number ( 15*,  (0SY(MN)/DSY(MN«2)  )  , *  1  •  J*0*0,2HI2) 

000650 

100 

CONTINUE 

000653 

CALL  PL0T(17.,((«YMTN/YSCALE)-,46)  ,-3) 

000660 

IF(NOTO*EO*NOT)  GO  TO  150 

000662 

NOTO  ■  NOTO«1 

000663 

GO  TO  1000 

000664 

150 

CALL  SYHBOL(0,0,0**.1«SH  ,90. ,5) 

000670 

CALL  PLOT (0,0,999) 

000673 

STOP 

000675 

END 

pnowku  length  including  i/o  duffers 
014200 


UNUSED  CDMPILER  SPACE 
014500 


N236565.  10/29/70, PUPOUF  MACE  70/10/26, 


18«45,1S,MZ365/  36634t JUOOvTPl tT60tCN60000tL30 

1S«45,18,00«P20. 

18«45,18,MAP(0N) 

18, 45,19, FUN(S) 

18.45.22,  CTIMC  001*316  SEC,  FllN  NOD  LEVEL  AOO 

18.45.22,  LOO. 

18. 45. 24,  CX  2.733  SEC,*  Px  3*31S  SEr. 

18. 45. 24,  NL  23500  WORDS 

18. 45. 24,  PLOTTING  STARTED 

18. 45. 42.  PLOTTING  SUCCESSFUL 
18, 45, 42. STOP 

18.45.43. COPYPLT(JUOO«36634tllx  90*2tBLACK) 

18. 45. 43.  CX  12.899  SEC.t  PX  8*752  SEC* 

18,45, 43*NL  4000  WORDS 

18,59,22*  MT51  ASSIGNED  -  CALCOMP 
19,00.25*MT51«  WPE  RECOVERED* 

19,01.35,COPYPLT  completed— tXPECT  PLOT  OUTPUT 
19,01.35,PLOT  0*602  HR, 

19,01,45*CP  20*832  SEC*.  PP  94*186  SEC* 

19,01.45,L1NES  1204 

19,01.45,CM  0*511  HWD-SEC,*  FL  4000  WoROS 


B15 


COMPUTER  PROGRAM  NO.  4 


00000? 

00000? 

000001 

000004 

000005 

oooon 

oooon 

0000?3 
000023 
000037 
000037 
00  ^054 
000056 
000060 
000105 
000105 
00011? 
000121 
000123 
000124 
000145 
000147 
000151 
000153 
000163 
000165 
00O166 
000200 
000210 
000210 
000212 
000213 
000220 
000221 
000223 
000233 
000234 
000241 
000241 
C00243 
000245 
000250 
000255 
000257 
000261 
000263 
000263 
000265 
000267 
0<00270 
000301 
000303 
000316 
000316 
000321 


PROGRAM  JUOD{INPUT*OUTPUT*PLOT.TAPE5«INPUT.TAPE6«OUTPUT*TAPE1«PLOT 
1) 

dimension  '40GR(5)  «CAL<8)  *047(140)  *HR(140)  *TIMIN(140)  *R(l40tS)  *7(14 
10*R) *0(140*8) *D5Y(140) •TX(140) 

CALL  PLOTS 
not  •  2 
NOTD  *  1 

1000  READ(5*1)  NA 

1  format (12) 

REAn(S*2)  NORfNOA 
2  FORMAT(/I3*I2) 

READ (5* 3)  OaY(1) ,HR(l) *TIMIN(1) tSN 
3  FORMAT(7X*3F5.0,7X*F5.0) 

READ (5* 4)  (CAL(J) * Jal *NOA) * (R ( 1 • J) t Jal *NOA) 

4  F0RMAT(/6X«8F8,n/6X*RF6,0) 

DO  5  K«2*N0R 

5  RFA0(5*6)  UAY(K) .HRCK) *TIH1N(K) * (R(K*J) »Jal*NOA) 

6  FORMAT(F4.0*FP,0*F4,0*PF8.2) 

TZERO  «  24.»0AV(1)  ♦  HR(1)  ♦  TIM!N(l)/60, 

DO  7  Jal*N0A 

7  7(1»J)  ■  0.0 
DO  9  Ka2*N0p 
no  A  Jal*NOA 

A  T{K*J)  ■  24.»0AY(K)  ♦  HR(K)  *  TIMIN (K) /6O.-T7ERO 
9  CONTINUE 
LL  ■  NA-1 
DO  10  JaltNA 
;0  D(1*J)  ■  0.0 

no  16  Ka2*N0R 

11  IF((R(K,NA),EQ.99999.).0R.(R(K*NA).E0. 00000*))  GO  TO  12 

D(K*NA)  a  (R(K*NA).R(l«NA))*rAL(NA) 

GO  TO  15 

1?  IF(K.EQ.NnR)  GO  TO  14 

NOP  ■  NOR.l 
DO  li  laKiNOR 
R(I*NA)  a  R(i*l,NA) 

T(I,NA)  a  TI^l.NA) 

DO  13  Lal.LL 
P ( I*L)  *  R( 1*1 *l  ) 

13  T(I*L)  a  T(I*lfL) 

GO  TO  11 

14  NOR  a  NOR-1 

15  NOGR(NA)  a  nor 

16  CONTINUE 

DO  17  Kal.Ll. 

17  N06R(K)  a  NOGR(NA) 

KNOA  a  NOA-1 

DO  50  JaltKNOA 
KOUNT  a  0 
NRR  a  NOGR(J) 

DO  50  Ka2*NBR 
KKEEP  ■  K 

IF((R(KtJ).CQ. 99999, ).0R*(R(K*J).EQ. 00000*))  60  TO  20 

IF(KOUNT*6T,0)  go  to  40 

0(K*J)  — (R(K*J)-R(i*J) )»CAL(J)*D(K*NA) 

60  TO  50 

20  IF(KKEEP,EQ,N0GR(J) )  GO  TO  45 
25  KOUNT  a  K0UNT*l 
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00324 

00320 

00327 

00340 

00341 

00343 

>00354 

'00376 

(00377 

100401 

100402 

100407 

00411 

S0412 
0414 
*00415 
00423 
00423 
00427 
100434 
00436 
[00437 
100441 
00442 
100447 
00447 
00453 
[00460 
(00463 
00465 
00474 
00474 
(00476 
00500 
|00516 
( 

[00516 

00523 

00537 

100541 

00541 

00542 

00543 

00544 

100551 

00555 

00561 

00563 

00565 

0057? 

00605 

00624 

00627 

00632 

00635 

00636 


30 

35 

40 

45 

50 


N06P(J)  ■  N06RIJ)-1 
NRR  ■  N06R(J) 

PO  35  KK«KKEEPfN8R 
P(KKfJ)  ■  R(Kk*1*J> 

T(KK*J)  ■  ^  aoo  .  UP  (D(KKEEP*J)  vFQ.OOOOO*) )  ©0  ^0  20 

fiO  TO  50 

NOGP(J)  •  N06R(J)-1 
igRR  ■  NOGR(J) 

CONTINUE 

YRAX  ■  0(lfl» 

no  200  J»1*N04 
(iIrR  ■  NOG*^(J> 

DO  200  I«2tNRR 

IF(n(I.J).BT.YMAX)  GO  TO  199 
60  TO  200 

199  yMAX  ■  0(I»J) 

200  CONTINUE 
YMIN  ■  D(l«l) 
no  202  J«l»NOA 
NRR  ■  N0G*^(J) 
no  202  I«2*NBR 
lF(n(I.J).LT,YMTN)  60  TO  201 

60  TO  202 
YMIN  ■  0(I»J) 

CONTINUE 

Y5CALE  ■  (YMAX-YMlNl/9. 

no  700  J«1«N06 

or  GOOD  W*0,NGS  .  .«3> 

NBR  •  NCG?- 

IlRlTE(6.7*:ii  «“<»J»/I“^’'^i!S‘5’lOX.7HTlME  ■  ♦FIO.?.IOX»15MOISPLACE 
FCRMATdH  tlOHREAOlNG  ■  fFlO,2f IOX.TMTIME  ■  triw.r 

IMENT  *  fF20«2) 

CONTINUE 
no  75  J«liNOA 
NN>N0GR(J)*1 
NNNbNN^I 
n(NNtJ)  ■  0, 

O(NNNfJ)  ■  YSCALE 
T(NNtJ)  ■  0. 

75  T(NNNfJ)  «  170*0  «  n  .O) 

CALL  SYMBOL(3,.l0*3*.l5t9HMPM  JO*:’®?’” 

CALL  NUHB€R(4,3,10*3*.l5*SNi0*#4Hr5.0) 

KN  ■  NOGR(l)  «1 

'?;a‘XMl';'MvlUsC*LE..UHTlHE  IN  HRS.-li .IS.G.O.O.T <KN. I .  .T .K 

1NN*1) f"0) 

CALL  PLOT (0* »0* **3) 

CALL  PLOT(15**0,**2) 

00  100  J«ltNOA 
MN«N00R( J) 

NOORJ«NOGR ( J) *2 


201 

202 


750 

751 
700 


BIT 


000640  DO  00  KeltNOGRJ 

000650  DSY(K)  s  DiKiJ) 

000651  RO  TX(K)  ■  T(K.J) 

000651  CALL  LINE (TX«OSYtMN»l tOtO) 

000657  CALL  NUMBtR ( lb. • (OSY (MN) /OSY (MN«2) ) » • 1 t J»0«0,2HI2) 

000666  100  CONTINUE 

000671  CALL  PL0T{17.» ( (♦YMTN/YSCALE)-,46) »-1) 

000676  IF(NOTU.EO.NOT)  60  TO  ISO 

000700  NOTO  ■  NOTO  ♦  I 

000701  60  TO  1000 

00070?  150  CALL  SYMBOL ( 0.0 1 0.  •  .1 *5H  tOO.tS) 

000706  CALL  PLOT (0«0f9O9) 

000711  STOP 

000713  END 

program  length  including  I/O  Buffers 
014220 

UNUSED  compiler  SPACE 
014400 


MUIIRRR.  10/27/70. PURDUE  MACE  70/10/26. 


19.09.19.MU11R/  36634, JUOO«TPlfT60fCM60000*L30 

19.09. 19, 00«P20. 

19.09. 19. MAP(ON) 

19.09.19.FUN(S) 

19.09.2P,CTIME  001.316  SEC.  FuN  MOO  LEVEL  60D 
19.09.P3.LGO. 

19.09.P8.CX  2.401  SEC.,  Px  2.532  SEc. 

l9.09.?fl.NL  23500  WORDS 

19. 09. 28. Plotting  started 

19. 10. 01, PLOTTING  SUCCESSFUL 

19,10.01, STOP 

19,10.02.COPYPLT( JUDO, 36634,1 lx  40,2tBLACK) 

19, 10. 02. CX  14,550  SEC,,  Px  i0.5T9  SEc. 

19, 10. 02. NL  4000  WORDS 
19,19.20,  MT52  ASSIGNED  -  CALCOMP 
19,19,25,MT52,  WPE  RECOVERED, 

19,22,?4.C0PYPLT  COMPLETED— tXPECT  PLOT  OUTPUT 
19, 22, 24, PLOT  0.294  MR, 

19. 22. 38.  CP  18.359  SEC.,  Pp  51.61*  SR. 

19. 22. 38.  LINES  2107 

19. 22. 38.  CM  0.466  MWO-SEC,*  FL  4000  WORDS 
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APPEHDIX  C 

This  drawing  replaces  Figure  3.27  in  Technicad  Report  No.  1 

(Programming  error  resulted  in  erroneous  displacement  values 
being  plotted  for  the  measured  data  on  the  original  figure.) 


Cl 


AXIAL  AXIAL  AXIAL 


(M3U!)  (M3U!) 

NOiSNaixa  N0ii3vtuN03  NOiSNaxxa  NoiiomNoa 
1VIXV  ivixv  1VIXV  ntix¥ 


C2 


b)  Left  Horizontal  MPBX  '  c)  Right  Horizontal  MPBX 


